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ABSTRACT 

Published B and V fluxes from nearby Type la supernovae are fitted to light-curve templates with 
4-6 adjustable parameters. Separately, B magnitudes from the same sample are fitted to a linear 
dependence on B — V color within a post-maximum time window prescribed by the CMAGIC method. 
These fits yield two independent SN magnitude estimates Rmax and Bbv- Their difference varies 
systematically with decline rate Amis in a form that is compatible with a bilinear but not a linear 
dependence; a nonlinear form likely describes the decline-rate dependence of Rmax itself. A Hubble fit 
to the average of Umax and Bbv requires a systematic correction for observed B — V color that can be 
described by a linear coefficient TZ = 2.59 ± 0.24, well below the coefficient Rb ~ 4.1 commonly used 
to characterize the effects of Milky Way dust. At 99.9% confidence the data reject a simple model 
in which no color correction is required for SNe that are clustered at the blue end of their observed 
color distribution. After systematic corrections are performed, Umax and Bbv exhibit mutual rms 
intrinsic variation equal to 0.074 ± 0.019 mag, of which at least an equal share likely belongs to 
Bbv- SN magnitudes measured using maximum-luminosity or CMAGIC methods show comparable rms 
deviations of order wO.Id mag from the Hubble line. The same fit also establishes a 95% confidence 
upper limit of 486 km s“^ on the rms peculiar velocity of nearby SNe relative to the Hubble flow. 
Subject headings: supernovae: general — cosmology: observations — distance scale 


1. INTRODUCTION 

Well before the Universe was found to be accel¬ 
erating llPerhi mtterel_aJ[IMSt iGajnaAdcheLaLl[19981: 
iSchinicit et al.lll998t iRiess et al.lL998t iPerlmutter et al.1 

Il999jl . Type la supernovae (SNe), the study of which un¬ 
derpinned that discovery, were thought to be standard 
candles. Even before the Type la subclass was identified, 
iKowall lll968^ s tudied the distrib ution of Type I SN lu¬ 
minosities, and iPskovskil l|1977ll observed that brighter 
SNe exhibit slower decline rates. The st a ndard ization 
of Type la SNe sharpened when iPhillinsI lll99,‘lll intro¬ 
duced a brighter-slower correction that was linear in the 
decline-rate parameter Amis (the change in magnitude 
from maximum light to 15 rest-f rame days the r eafter) . 
That c orrection was r efined by iHamuv et al.l l|1996a|l 
and by iPhillins et al.1 (h encefo r th Phil99), who 

added a quadratic term. Trind l[l998ll . iTriop fc Branchl 
ifT^ . and iParodi et al.l (1200611 fit a correction that 
was linear both in Amis and in a second parameter, 
B — V color. In these two-parameter studies the ob¬ 
served brighter-bluer correlatio n was compensated by a 
fit c olor coefficient TZ » 2.5 . IPerlmutter et al.1 lll997ll 
and IColdha.ber. Groom et al.l 1)200111 thenceforth Perl97 
and GoldOl) stretched the time axis of a fixed template to 
approximate different SN light curves; that linear stretch 
factor substituted for Amis in their alternative brighter- 
slower correction. 

Given an adequate training set - here an ensemble of 
SNe in the nearby Hubble flow whose deviations A^ from 
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a fiducial absolute SN magnitude are deduced from their 
redshifts - the problem of estimating the deviation A of 
a single SN, using a set {xj} of its measured parame¬ 
ters, is amenable to general solution. In principle the Xj 
can include any parameter that may add to knowledge 
of A: fluxes in various bands at various phases, spectral 
line amplitudes, etc. One example of such a solution is 
given in Appendix^ Present progress along these more 
general lines is exemplified by t he MLGS2k2 package and 
its pr edecessors described bv iRless. Press, fc Kirshn^ 
lIT^ IPiess et al.l (l199sll . and l.lhal (1299211 . ~ 

those authors found, multivariate approaches to SN stan¬ 
dardization are severely limited by low training statis¬ 
tics. For example, the MLCS2k2 training set did not ade¬ 
quately constrain the color-correction coefficient i?y, for 
which an external value was imposed. 

This paper aims to extend understanding of the sys¬ 
tematic corrections that help standardize Type la super¬ 
nova luminosities, and of the intrinsic variations which 
blur that standard. Large resources now are devoted 
to projects that use standardized SNe to measure with 
exquisite precision the parameters determining the Uni¬ 
verse’s recent acceleration, and even larger efforts will 
be needed to track the evolution of its equation of state, 
revealing the sources of that acceleration. Even small ad¬ 
vances in such understanding can have big implications 
for the realization of these projects. 

In this paper, we make use of only a few properties of 
nearby SNe - luminosity, color and decline rate - that 
are simple to characterize. In that sense our approach 
to SN standardization is traditional, in the vein of work 
cited in the first paragraph. Nevertheless, we do bring 
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fresh tools to this study. By taking the difference be¬ 
tween the two simplest measures of SN magnitude - lu¬ 
minosity at peak, and luminosity at hxed color using the 
CMAGIC method llWang et al.ll200,‘lt henceforth WangOS) 
- we impose a new, more precise constraint on the (non¬ 
linear) form of their decline-rate dependence, and we set 
a robust lower bound on the intrinsic variation suffered 
by at least one of these measures. By applying a newly 
extended light-curve fitting method and standard error 
propagation consistently to each SN in the sample, we 
obtain a Hubble ht that sets a significant upper bound 
on the rms peculiar velocity of SNe relative to the Hubble 
flow. 

In Section |21 our light-curve (lc) and CMAGIC fitting 
methods are introduced. Their outputs are shown to 
be independent, and the corrections made to them for 
extinction and time dilation are described. Section |3 
identifies the photometric data sample that is processed 
by these methods, and the subsample of SNe to which 
Hubble hts are applied. Section 0 reveals the new in¬ 
formation gained by comparing the LC and CMAGIC ht 
magnitudes. In Section [3 the result of a Hubble ht to 
the average of these outputs is presented; this is the h- 
nal step in a global ht to the sample. Section El uses the 
global ht parameters to quantify the deviations of the 
LC and CMAGIC outputs from the Hubble line. Section 
0 explores properties of the subset of SNe assumed to 
be largely free of extinction by host-galactic dust. The 
paper concludes with a summary and discussion. 

2. METHODS FOR FITS TO PHOTOMETRIC 
DATA 

We used two complementary methods to estimate SN 
magnitudes. The hrst, a light-curve method described 
in EH is based on the ht peak magnitudes in the B 
and V bands, and on the ht diherence Amis between 
the peak B magnitude and its value 15 rest-frame days 
later. These three key numbers and their uncertainties 
are supplied by our own LC hts to published measured 
huxes in these bands. The hts were performed by a newly 
and substantially extended version of the algorithm of 
Perl97 and GoldOl, in which the time axes of hxed tem¬ 
plates were linearly stretched to best represent the hux 
data. After color and decline-rate corrections, the SN 
magnitude yielded by our LC method is called Rmax- 

Our second meth od for estimating SN magnitudes, 
summarized in E2 is based on the output of the CMAGIC 
algorithm introduced by WangOS. This algorithm ex¬ 
ploits the nearly linear dependence of B magnitude upon 
B — V color over a period extending, for typical SNe, 
from «6 to f«27 rest-frame days after maximum light. A 
least-squares linear ht to published magnitudes and col¬ 
ors measured within that range yields an intercept Bbv 
at hxed color B — V = 0.6. The above mentioned LC ht in 
the B band is also used to supply Ami 5 and the phase at 
maximum. After color and decline-rate corrections, the 
SN magnitude yielded by our CMAGIC method is called 

Bbv- 

By no means are either of these methods restricted 
to the B and V bands. Much study has been devoted to 
their extension to the R, /, and U bands (WangOS; Wang, 
L. et al. 2005, in preparation, henceforth Wang05a; Gold- 
haber, G. et al. 2005, in preparation). However, since this 
paper relies on comparing the results of these two meth¬ 


ods, sufficient data must exist to complete both types 
of ht for each SN. Were we additionally to require suc¬ 
cessful LC and CMAGIC hts to the same SNe in the R, I, 
or U band, our sample size at best would halve. Data 
in these additional bands therefore remain outside this 
paper’s scope. 

2.1. Light Curve Fits 

This study requires use of a light-curve ht algorithm 
only for interpolating or extrapolating measured huxes 
to compute Hmax, Vniax, and Ami 5 . Here no other use 
is made of the LC ht parameters, and no physical sig- 
nihcance need be ascribed to them. For completeness, 
neverthess, we describe briehy the parametrization that 
we employed. More information on this algorithm is sup¬ 
plied in Wang05a. 

The set of published hux measurements (for example 
in the B band) for a single SN was ht to the function 

/(t; To, to) So, t-q, eo, td) = ro[T^(t*(< — to;so,'ro))-|- 

eoTBit*{t-to-td;so,To))] ( 1 ) 

where t is the rest-frame phase of the SN hux measure¬ 
ment; To is an overall normalization; TB{tr) is a hdu- 
cial template for B band flux vs. rest-frame phase G 
that peaks at G = 0; and to is the phase at which the 
hrst (main) term reaches its peak. In the hrst term, the 
“stretched” phase t*, discussed further in Appendix IHl is 
an odd function of t — to with parameters so and tq. The 
“stretch” parameter so is equal to dt/dt* at t = to; the 
“duration” parameter tq is such that dt/dt* is damped to 
unity when |t — to| 3> tq. Relative to the main term, the 
second (delayed) term is different only in that its phase 
at peak amplitude is delayed by td and its amplitude is 
reduced by the factor eg- The special case {tq = oo (no 
damping); eo = 0 (no delayed term)} was studied by 
Perl97 and GoldOl; under those conditions their param¬ 
eter s is equivalent to so. In all hts here, the parameter 
To is allowed to vary along with Tq, to, and sq; when eo is 
nonzero, as is permitted (Appendix 0]) for a minority of 
SNe, eo and td vary as well. Within a band, usually all 
available data were htted, including those in the “tail” 
region t* > 30 dy. Here we are not concerned with re¬ 
sults that are specihc to the tail region, or the transition 
thereto, because we emphasize methods that may prof¬ 
itably be applied to faint SNe that (unlike the nearby 
sample) are sensitive to cosmological acceleration. 

For neighboring bands the LC ht procedure is the same 
(except for choice of template). Smoothed versions of the 
B template d escribed by GoldOl and of the V template 
described by iKnop et'ffl ll2003ll are used. Parameters 
that are fit to data in different bands are nearly indepen¬ 
dent, with only two exceptions: for all bands in which a 
given SN is fitted, eo either is left free or is set to zero; 
and the epochs of peak hux are required to align within 
certain loose windows. To satisfy the latter conditions, 
all bands were fit at once. For decline-rate corrections 
we did not use Amj^ in the V band as a substitute for, 
or in combination with, the B band value: though a sat¬ 
isfactory correlation of Am}); with Ami 5 was seen, and 
reasonable Hubble fits using Amj}; alone or in combina¬ 
tion with Amis were produced, tighter Hubble residuals 
were obtained by using Amis alone. Neither did we sub¬ 
stitute the fit parameter I/sq for Amis, because Amis 
is correlated also with the companion parameter tq. 
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For various choices of the stretch function t* {t) in equa¬ 
tion o, our LC fit results for i?max and t4iax are sta¬ 
ble. Measuring Arn-is is more delicate: because the time 
derivative of the flux at day 15 can be large, a small error 
in the definition of the date of maximum flux can trans¬ 
late into a large error in Amis. Especially near those 
times, both t*{t) and the template parametrization must 
be smoothly differentiable. Accordingly we employed a 
particular choice of t*{t) for which details are supplied 
in Appendix 0 

In summary, we have applied to the full SN sample 
a uniform mathematical approach for determining their 
peak magnitudes and decline rates. The approach is self- 
consistent, with covariance matrices used to propagate 
the raw flux errors through to the final uncertainties. As 
a check of these methods, we subjected accurate graphs 
of each of our LC fits to detailed inspection which con¬ 
firmed that, within quoted uncertainties, the measured 
fluxes are represented satisfactorily by these fits. 

As an alternative approach, we could have attempted 
to combine and use the published results of different 
light-curve fits, employing various methodologies, that 
had been performed on these SNe. The dangers in 
such an attempt become evident when one compares 
our values of parameters such as Amis and peak mag- 
ni tude with th o se fou nd e.g. in the useful compilation 
of IReindl et al.l II2005D . Significant differences are found, 
particularly in values of Amis. For example, the com¬ 
piled values of Amis are largely confined above a lower 
limit of « 0.84-0.87 mag, with several exactly equal to 
0.87. Other SNe are clustered at Amis = 1-13 and 1.69. 
Such artifacts, which are not found in our own Amis 
distributions, appear to have arisen from limitations in 
the number of well-observed SN light curves that were 
available for use as discrete templates by the published 
fits. 

2.2. Fits to B Magnitude vs. B — V Color 

The CMAGIC method described by Wang03 begins with 
pairs of fluxes measured in the B and V bands at the 
same epoch U. These flux pairs are converted to raw 
magnitude pairs {Bi,Vi} and entered on a plot of Bi 
vs. {Bi — Vi). In parallel, the above-described LC fit is 
used to determine the decline-rate parameter Amis and 
the rest-frame phase tpeak of peak B flux. To those pairs 
{Bi,Vi} which lie within a fiducial region tf, — tpeak < 
ti < te — tpeak, a weighted least-squares linear fit is per¬ 
formed in which the errors on Bi and Vi are taken to be 
independent. At least three pairs must participate in the 
fit, providing a minimum of one constraint. The fit out¬ 
puts are B'^y, the intercept of the best-fit line with the 
fixed color Bi — Vi = 0.6 mag; the slope (}bv of that line; 
and their err ors. (Our BVy is e quivalent to Bbvo.6 of 
Wang03 and iConlev et all (1200,511 (henceforth Conl05).) 
The boundaries {tf, — tpeak) te — tpeak} of the CMAGIC lin¬ 
ear region are {(7 mag)/Amis, (30 mag)/Amis} dy for 
normal SNe, and {(10 mag)/Amis, (29 mag)/Ami 5 } dy 
for spectroscopically SN1991T-like SNe. 

As noted by Wang03, the average slope Pbv is very 
nearly equal to 2, with an rms intrinsic scatter (apart 
from photometric error) that is less than 0.1. If I3bv ei¬ 
ther were assumed or were measured always to have the 
same value, as is nearly the case, our choice of intercept 
color Bi — Vi = 0.6 would be immaterial. Here this par¬ 


ticular value is chosen for compatibility with earlier work 
(Wang03; Conl05), and for its consistency with the inter¬ 
cept color that minimizes both the measured rms Hubble 
residual of Bbv and the | covariance | between the inter¬ 
cept and [3bv- 

2.3. Photometric Errors 

Sources of photometric SN data are diverse; reported 
errors on data points are not always reliable, and their 
covariances normally are unavailable. Additional system¬ 
atic errors arise from observations on different telescopes 
and from uncertainties in the adopted light curve mod¬ 
els. To gain confidence in the errors on LC and CMAGIC 
fit outputs, we examined the distributions of confidence 
levels (CLs) for each type of fit. Ideally both distributions 
should be uniform between 0 and 1. However, typically 
the LC fits had smaller CLs than expected, likely because, 
even with 4-6 adjustable parameters, the templates were 
not flexible enough to fully track the detailed evolution of 
well-measured flux data. Conversely, the CLs for CMAGIC 
fits typically were larger than expected, perhaps because 
published photometric errors in the CMAGIC linear re¬ 
gion are slightly conservative. To address this issue for 
each type of fit, we mapped the observed CL distribution 
onto a uniform one while maintaining each SN’s rank in 
the distribution. This was done by scaling all fit errors 
that were output for a particular SN by (XOTig/Xunif)^^^’ 
where Xorig ^^e original x^ statistic and xSnif i® 
value which allows that SN to maintain its rank in the 
uniform CL distribution. For the LC and CMAGIC fits, the 
median error-scaling factors were 1.57 and 0.78, respec¬ 
tively. 

2.4. Extinction o/ Hmax and Bbv 

As a first step, standard corrections for the effects of 
Milky Way (mw) dust were applied to all of the LC and 
CMAGIC fit outputs that require it. This step will be 
discussed in m In this section, on the other hand, we 
are concerned with extinction by intervening dust outside 
the MW. An experimentally related issue is the effect 
upon SN luminosity of variation in intrinsic SN color. 

If intervening dust outside the MW causes true B — V to 
differ from raw {i.e. observed) (B — V)'^™ by 5{B — V) = 
(B — V) — (i? —E)™ , then, according to the law of 
ICardelli et al.1 II1989I1 (henceforth Card89), all true B 
magnitudes differ from raw B magnitudes by SB = 
^_^raw _ Rb 5{B — V), where Rb is approximately 
constant and of size ~ 4.1 if the dust is parametrized 
as is common for MW dust. For maximum-luminosity 
analysis, the SN color is measured here by the differ¬ 
ence E = between the B magnitude at 

7?max and the V magnitude at t4iax, as determined, re¬ 
spectively, by the nearly independent LC fits in the B 
and V bands. More generally, using the Card89 law to 
describe the effects on raw SN luminosity both of in¬ 
trinsic S N color variation and of intervening dust, and, 
following ivan den BerghI l|1995ll . denoting by 7Z the color- 
correction coefficient parametrizing those combined ef¬ 
fects, color-corrected Hmax is related to color-uncorrected 
SS,Tx by 

Smax = S(nTx ” ^ZE . (2) 

(Beginning in 21 we shall use the symbol Hmax to de¬ 
note the maximum B magnitude after full corrections. 
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which depend both on observed color and on decline rate, 
are applied. In the present discussion, nevertheless, for 
simplicity we do not include a decline-rate correction to 

For CMAGIC analysis, however, equation © does not 
hold: the change due to intervening dust in the CMAGIC 
output Bbvi defined as the B magnitude at fixed color 
B — V, will be smaller. Because it must continue to be 
measured at the same fixed color, the true value of Bbv 
should differ from the raw value B'^y by only that por¬ 
tion of 5B which is orthogonal to S{B — V): 

= {Rb- Pbv)S{B-V) . (3) 

Again taking into account the effects both of intrinsic SN 
color variation and of intervening dust, the corresponding 
effective color-correction coefficient for Bbv is 7?. — (5bv ■ 

In raw form, before correction for color or decline rate, 
B'^y uses only the B and V fluxes that are measured in 
the CMAGIC linear region, beginning at least several days 
after maximum light. Nevertheless, to define that region 
and to measure the decline rate, in parallel one performs 
a light-curve fit that does also use the fluxes measured 
near maximum, normally in the B band. To measure 
and correct for the SN color, in principle a second LC 
fit in the V band would allow the maximum-luminosity 
color E to be calculated. However, to preserve great¬ 
est orthogonality to maximum-luminosity analysis, the 
CMAGIC method substitutes a complementary measure 
of SN color that is unaffected by the details and uncer¬ 
tainties of the V light curve near maximum. This is done 
by exploiting the unique linear combination 

Draw _ Draw 

£ — ^max ^BV 

Pbv 

of and that permits S£ = 5{B — V) in the 
example to which equation m pertains. For CMAGIC 
analysis, £ is not the only possible color measure, but it 
is a choice that has become conventional and is adopted 
throughout this paper. 

In analogy to equation the arguments following 
equation © require 

Bbv = - ( 7 ^ - Pbv)£ ■ ( 5 ) 

Using equation 0 to evaluate £, equation 0 becomes 

BBv = Bi:Z.-n£ . ( 6 ) 

Comparison of eqs. and 0 emphasizes that Bbv 
and Hniax are independent only in that they use different 
measures of SN color to correct When the CMAGIC 

color £ closely tracks the maximum-luminosity color E, 
as is usually the case, Bbv and Hmax maintain the same 
sensitivity to uncertainties in TZ. 

2.5. Covariance of By^ax and Bbv 

Given the close relationship between Hmax and Bbv 
revealed by eqs. m and it is natural to examine 
their mutual correlation. Using a Monte Carlo technique, 
Conl05 studied the covariance due to photometric error 
between B'^y and after correcting these quantities 

only for decline-rate dependence. In that study, 
was determined by an LC fit of the type used by Perl97 


and GoldOl, as discussed in A mean Pearson cor¬ 

relation coefficient of «0.15 was determined. Since a co- 
variance at that level would not have had a significant 
effect, correlations due to photometric error between en¬ 
tirely uncorrected B’^^y and were not taken into 

account in the error propagation performed here. On 
the other hand, correlations due to photometric error 
in corrections to Bbv and By^yx are more significant: for 
example, after extinction corrections, both quantities de¬ 
pend on B^^^. Such correlations are included straight¬ 
forwardly in our error propagation, as are correlations 
due to peculiar velocity and extinction by MW dust. 

On a related note, it is of possible academic interest to 
identify the conditions, if any, that would cause corrected 
Bbv and Hmax to become fully correlated. We find that 
when Pbv is fixed, and stretch and time-of-maximum are 
determined for both B and V bands by a H-template 
fit, the maximum-luminosity and CMAGIC methods are 
fully correlated if “magic templates” are used and if V 
observations outside the CMAGIC linear region are dis¬ 
carded; and if the number of B observations is minimal, 
or if the B points used by CMAGIC are taken from the 
B template fit. (“Magic templates” uib and my satisfy 
diriB/<i{mB — rnv) = Pbv in the CMAGIC linear region.) 
None of these conditions are met by the analysis reported 
here. Appendix IgI demonstrates that Bbv and Hmax in¬ 
deed are fully correlated for a hypothetical case that does 
satisfy these requirements. 

2.6. K Corrections 

The Hubble expansion not only dilates the apparent 
SN evolution, but also causes its light to be bluer at 
rest-frame emission than at detection. To translate pub¬ 
lished photometric data from observed to standard rest- 
frame passbands, we employed a particular AT-correction 
method, discussed more fully in Wang05a, that is de¬ 
scribed only briefly below. 

For the actual range of redshifts 0.003 < z < 0.079 
of the SNe in our fitted sample, light detected at the 
peak of the V (B) passband was emitted at a rest-frame 
wavelength that is blueshifted by up to 35% (60%) of the 
wavelength interval by which the peak of the adjacent B 
(U) passband is separated. Nevertheless, in part because 
the light curve is fit nearly independently in each pass- 
band, the functional forms used by both our LC (EU 
and CMAGIC 1 ^2.211 fits are flexible enough to permit both 
types of fit to be performed to the observed data before 
K corrections for these wavelength shifts were applied; a 
posterior K correction needed to be applied only to the 
five fit outputs B^yff^yy, Knaxi , and (3bv- We 

index these outputs by k {l<k<5). 

The framew ork for our K correc tions begins with that 
established bv iHamuv et al.1 lll993ll . who directly applied 
the spectrophotometric data from three library SNe to 
the task of AT-correcting a time series of observed B and 
V magnitudes belonging to other SNe that lacked spec¬ 
tral data. At present we benefit from a library of >30 
SNe for which spectroscopic as well as good photomet¬ 
ric data are available. For the B and V bands we con¬ 
structed a set of SN light curve pairs, with each pair 
characterized by a different raw color Ei. We then used 
the spectroscopic data to modify each pair so that it cor¬ 
responded to each of a set of hypothetical redshifts Zj. 
After applying standard LC and CMAGIC fit procedures to 
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TABLE 1 

Number of SNe passing cumulative cuts that yield the 

BASE, PLOTTED, AND FITTED SAMPLES. 


Condition satisfied (cumulative) 

No. of 
SNe 

retained 

Sample 

name 

No. of SNe initially considered 

110 


Enough data to fit Bmax and Bbv 

84 


z < 0.1 

82 


Amis < 1.7 

77 


Bmax — VmEix 0.5 

75 


1.5 < dsv < 2.5 

71 


Not peculiar (SN2002cx) 

70 

base 

(B-w )(VA-w)+<^2 < 0.25 

61 

plotted 

Each SN adds < 8 to global fit 

56 

fitted 



Fig. 1.- SNe per 0.025 vs. S = + 

The arrow shows the position of the cut S < 0.25. 


the entire collection of light curve pairs, we recorded the 
differences {Kk{Ei, zj)} between fit output k at z = Zj 
and the same output at z = 0. For a given j and k, the 
corrections {Kk{Ei, Zj)} were fit to a linear function of 
the Ei, with parameters fj,j and In turn, the {/tj } 
and each were fit to cubic polynomial functions of 

the Zj. When applied to a real SN, the K correction to 
the k^^ fit output was calculated by evaluating the poly¬ 
nomial that was fit to the {Kk{Ei, Zj)} at that SN’s raw 
color E and redshift z. 

We found this scheme attractive because its imple¬ 
mentation is straightforward and because the K correc¬ 
tions are smooth functions of z and, like Bi itself in the 
CMAGIC region, they are linear functions of color. Exam¬ 
ples of results from this scheme are provided in Appendix 

ini 

3. FITS TO PHOTOMETRIC DATA 

For inclusion in this analysis, the available sample of 
nearby SNe was subjected to a set of standard cuts that 
are listed in Tabled First, as noted in m enough data 
were needed so that both a light-curve and CMAGIC fit 
could be made to each SN. Primarily this meant that at 
least three pairs {Bi,Vi} of magnitudes were required 
to have been measured in the CMAGIC linear region. 
These requirements eliminated one quarter of the can¬ 
didates. Secondly, to avoid sensitivity to cosmological 
acceleration and to bound the range of redshift over 
which K corrections are performed (fOl), we required 
z < 0.1. Thirdly, we consider highly subluminous and 
fast-declining “SN1991bg-like” SNe to be a recognizably 
distinct subset. If lumped together with the main sam¬ 
ple, they would only muddy its characteristics. As this 
«6% subset is too small to analyze separately, we sim¬ 
ply cut it out by requiring Amis <1-7 mag. Fourthly, a 
standard loose cut A < 0.5 mag on raw B — V color elim¬ 
inated only two severely reddened candidates. As noted 
in a the CMAGIC fit slope Pbv is clustered near 2 with 
an intrinsic scatter <0.1. Our fifth cut, 1.5 < (5bv < 2.5, 
eliminated four SNe whose CMAGIC slopes are patholog¬ 
ical and/or poorly measured. Finally we discarded a 
single candidate, SN2000cx, whose properties are widely 
considered to be highly atypical of SNe la. (Other candi¬ 
dates that also might be considered to be peculiar failed 
earlier cuts.) These standard cuts left a “base sample” 
of 70 SNe. 


To the base sample we applied a “well-measured” cut 
which requires explanation. Normally, when several mea¬ 
surements of a physical quantity are combined and a 
proper error analysis is applied, nothing is lost by in¬ 
cluding measurements that are far less precise than those 
which dominate the result. One of our main interests 
here, however, is the intrinsic variation of absolute SN 
magnitudes. In the simplest picture, this is determined 
by subtracting the calculated variance due to photomet¬ 
ric error from the observed total variance. If unusually 
poorly measured SNe are included, both variances in¬ 
flate, and the precision of their difference suffers. As will 
be seen in SH our main information on intrinsic SN vari¬ 
ation comes from comparing the maximum-luminosity 
B — V color E to the CMAGIC color £. As is discussed in 
Appendix H the photometric error on AE = £ — E is 
tracked roughly by the simple quantity 

S = , (7) 

where the ct’s are raw errors output by the LC and 
CMAGIC fit procedures described in io and iO To 
identify a stable sample for further analysis, we chose to 
cut on S, as opposed to the precise error on AE, because 
the latter depends on fit parameters that were expected 
to evolve as the analysis matured. Figure ^ exhibits the 
distribution in E of the 70 SNe in the base sample. Re¬ 
maining after the cut S < 0.25 are the 61 SNe which 
comprise the “plotted sample”. 

For completeness. Tabled lists a final cut that will be 
justified in ^5.31 Unlike all cuts discussed so far, the 
final cut takes into account the residuals from fits that 
demand consistency between maximum-luminosity color 
E and CMAGIC color £ (0. and from Hubble fits to the 
average of Hmax and Bbv (0- It excludes as outliers 
five of the 61 SNe in the plotted sample, yielding a “fitted 
sample” of 56 SNe. These outliers are excluded from the 
AE = 0 and Hubble fits, but they continue to be included 
in the figures, identified by open-circle symbols. 

Table 121 presents the results of our LC and CMAGIC fits 
to the plotted sample. Tabulated are the SN name; the 
symbol if the host galaxy is E, E/SO, or SO; the num¬ 
ber of adjustable parameters in the LC fit to one band 
(see O; the raw peak luminosities in the B and V 
bands; the decline rate Amis; the raw CMAGIC luminos¬ 
ity B'^y and slope I3bv', the redshift z in the CMB frame; 
the excess B — V color E’^'^(H —U) due to Milky Way 
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TABLE 2 

Results of light-curve and cmagic fits. 


SN 

E or 
SO 
host 

No. of 
params 
in LC fit 

Draw 

•^max 

(mag) 

■p^raw 

^max 

(mag) 

Amis 

(mag) 

Draw 

fBV 

(mag) 

0BV 

•^cmb 

£MW 

(B-V) 

(mag) 

Ref 

1990N 


4 

12.537(031) 

12.626(011) 

1.134(023) 

13.831(026) 

2.123(113) 

0.00447 

0.026 

(1) 

19900 


4 

16.223(056) 

16.248(035) 

0.937(059) 

17.575(021) 

2.118(059) 

0.03060 

0.093 

(2) 

1990af 

★ 

4 

17.767(018) 

17.728(012) 

1.594(030) 

18.965(052) 

2.038(190) 

0.04992 

0.035 

(2) 

1991T 


6 

11.412(054)'" 

11.455(009) 

1.183(084)'" 

12.415(082) 

2.029(236) 

0.00693 

0.022 

(3) 

1991U 


4 

16.487(102) 

16.577(053) 

1.172(044) 

17.646(029) 

1.961(079) 

0.03096 

0.062 

(2) 

1991ag 


6 

14.452(060) 

14.419(073) 

0.896(033) 

15.807(050) 

1.982(116) 

0.01388 

0.062 

(2) 

1992A 

★ 

4 

12.609(018) 

12.489(018) 

1.349(017) 

13.654(040) 

2.006(133) 

0.00594 

0.017 

(4) 

1992ag 


6 

16.250(058) 

16.152(033) 

1.080(066) 

17.226(135) 

1.842(257) 

0.02734 

0.097 

(2) 

1992al 


4 

14.479(017) 

14.531(014) 

1.217(017) 

15.923(012) 

2.179(032) 

0.01350 

0.034 

(2) 

1992bc 


6 

15.107(039) 

15.142(019) 

0.750(024) 

16.816(035) 

2.074(116) 

0.01978 

0.022 

(2) 

1992bg 


4 

16.753(046) 

16.770(026) 

1.160(039) 

18.097(031) 

2.189(092) 

0.03648 

0.185 

(2) 

1992bh 


4 

17.614(031) 

17.546(015) 

0.972(078) 

18.730(033) 

2.112(092) 

0.04530 

0.022 

(2) 

1992bk 

★ 

4 

18.259(068) 

18.214(049) 

1.631(043) 

19.377(028) 

1.720(122) 

0.05885 

0.015 

(2) 

1992bl 


4 

17.409(055) 

17.374(040) 

1.502(045) 

18.661(060) 

2.276(183) 

0.04223 

0.011 

(2) 

1992bo 

★ 

4 

15.748(021) 

15.754(019) 

1.631(019) 

16.925(027) 

1.838(098) 

0.01811 

0.027 

(2) 

1992bp 


4 

18.258(029) 

18.308(022) 

1.354(031) 

19.605(056) 

1.721(153) 

0.07858 

0.069 

(2) 

1993B 


4 

18.437(062) 

18.431(038) 

1.321(047) 

19.610(055) 

1.867(194) 

0.07009 

0.079 

(2) 

1993H 


4 

16.709(023) 

16.528(018) 

1.566(026) 

17.403(033) 

2.039(090) 

0.02507 

0.060 

(2) 

19930 

★ 

4 

17.615(025) 

17.672(021) 

1.278(027) 

18.966(035) 

2.149(187) 

0.05293 

0.053 

(2) 

1993ag 

★ 

4 

17.803(028) 

17.722(023) 

1.283(039) 

18.861(049) 

1.983(243) 

0.05004 

0.112 

(2) 

1994D 

★ 

4 

11.762(019)'" 

11.855(019)'" 

1.555(027)"" 

13.223(048) 

2.089(125) 

0.00261 

0.022 

(5) 

1994M 

★ 

4 

16.324(034) 

16.265(021) 

1.409(033) 

17.500(111) 

1.993(454) 

0.02431 

0.024 

(6) 

1994Q 

★ 

4 

16.518(108) 

16.428(060) 

0.961(061) 

17.599(041) 

2.033(200) 

0.02897 

0.017 

(6) 

1994S 


4 

14.808(029) 

14.807(041) 

0.869(099) 

16.208(026) 

1.989(067) 

0.01616 

0.021 

(6) 

1994ae 


4 

13.117(055)'" 

12.985(032) 

1.072(052)"" 

14.402(028) 

2.050(122) 

0.00539 

0.031 

(6) 

1995D 

★ 

4 

13.230(016) 

13.264(014) 

0.977(015) 

14.610(030) 

2.178(117) 

0.00654 

0.058 

(6) 

1995ac 


6 

17.056(032)'" 

17.094(032)'" 

0.852(070)"" 

18.444(026) 

2.050(127) 

0.04882 

0.042 

(6) 

1995ak 


4 

16.056(046) 

16.056(047) 

1.486(054) 

17.184(038) 

1.929(174) 

0.02300 

0.043 

(6) 

1995al 


4 

13.319(017) 

13.187(018) 

0.964(017) 

14.318(038) 

2.027(116) 

0.00587 

0.014 

(6) 

1995bd 


6 

15.220(016) 

14.910(016) 

0.776(035) 

16.172(062) 

2.163(093) 

0.01599 

0.495 

(6) 

1996C 


4 

16.617(077)'" 

16.528(041)'" 

0.940(038)"" 

17.856(040)'" 

2.182(115)'" 

0.03007 

0.014 

(6) 

1996X 

★ 

6 

12.978(018)'" 

13.007(018)'" 

1.307(021)"" 

14.271(017) 

2.037(042) 

0.00800 

0.069 

(6,7) 

1996bl 


4 

16.677(029) 

16.627(021) 

0.851(038) 

17.938(024) 

1.928(068) 

0.03485 

0.105 

(6) 

1996bo 


4 

15.838(023) 

15.515(012) 

1.171(030) 

16.346(021) 

2.096(091) 

0.01632 

0.078 

(6) 

1996z 


4 

14.409(091) 

13.992(098) 

1.060(073) 

14.901(022) 

2.487(387) 

0.00796 

0.063 

(6) 

1997E 

★ 

4 

15.108(019) 

15.055(013) 

1.394(028) 

16.233(014) 

1.836(062) 

0.01332 

0.124 

(8) 

1997bp 


4 

13.891(026)'" 

13.761(026)'" 

1.288(109)"" 

14.651(065)'" 

1.831(121)'" 

0.00944 

0.044 

(8) 

1997bq 


4 

14.442(025) 

14.342(024) 

1.160(020) 

15.276(018) 

1.994(049) 

0.00960 

0.024 

(8) 

1997br 


6 

13.596(028)'" 

13.418(028)'" 

1.139(034)"" 

14.026(072) 

2.313(148) 

0.00630 

0.113 

(9) 

1997CW 


6 

15.863(067) 

15.586(081) 

0.997(034) 

16.488(030) 

1.823(083) 

0.01595 

0.073 

(8) 

1998ab 


6 

16.137(076)'" 

16.113(076)'" 

1.098(035)'" 

17.200(024) 

2.247(106) 

0.02784 

0.017 

(8) 

1998bu 


6 

12.130(015) 

11.774(010) 

0.977(024) 

12.842(039) 

2.015(088) 

0.00416 

0.025 

(10) 

1998dm 


4 

14.761(042) 

14.551(029) 

0.986(032) 

15.392(036) 

2.146(115) 

0.00556 

0.044 

(8) 

1998ec 


4 

16.309(059) 

16.134(056) 

1.062(031) 

17.088(027) 

2.051(079) 

0.02012 

0.085 

(8) 

1998es 


6 

13.836(013) 

13.760(008) 

0.742(015) 

15.099(031) 

1.654(177) 

0.00957 

0.032 

(8) 

1999aa 


6 

14.741(014) 

14.722(009) 

0.810(019) 

16.165(008) 

1.801(029) 

0.01525 

0.040 

(11) 

1999aw 


6 

16.694(028)'" 

16.745(019)'" 

0.787(026)'" 

18.245(015) 

1.736(071) 

0.03924 

0.032 

(12) 

1999dk 


4 

14.817(027) 

14.783(014) 

1.002(029) 

15.880(025) 

2.148(076) 

0.01395 

0.054 

(13) 

1999dq 


6 

14.405(016) 

14.342(009) 

0.953(028) 

15.475(055) 

2.074(106) 

0.01308 

0.110 

(8) 

1999ee 


6 

14.853(007) 

14.566(007) 

0.904(008) 

15.650(014) 

2.081(032) 

0.01055 

0.020 

(14) 

2000E 


6 

12.784(064)'" 

12.757(064)'" 

1.070(115)'" 

13.884(010) 

1.752(015) 

0.00422 

0.364 

(15) 

2000ce 


4 

16.961(061) 

16.518(030) 

0.961(032) 

17.152(032) 

2.183(059) 

0.01650 

0.057 

(13) 

2000cf 


4 

17.087(059) 

17.103(018) 

1.316(056) 

18.327(034) 

2.068(111) 

0.03646 

0.032 

(8) 

2000cn 


4 

16.601(009) 

16.399(010) 

1.582(027) 

17.322(036) 

1.972(389) 

0.02321 

0.057 

(8) 

2000dk 

★ 

4 

15.354(027)'" 

15.330(012) 

1.440(042)'" 

16.424(030) 

1.971(157) 

0.01645 

0.070 

(8) 

2001V 


6 

14.606(064)'" 

14.546(064) 

0.730(081)'" 

15.766(086) 

1.858(286) 

0.01604 

0.020 

(16) 

2001ba 


4 

16.322(025) 

16.342(014) 

0.911(024) 

17.761(059) 

2.367(255) 

0.03053 

0.064 

(17) 

2001el 


4 

12.793(027)"" 

12.691(027)'" 

1.161(022)'" 

13.621(008) 

2.092(028) 

0.00365 

0.014 

(18) 

2002bo 


4 

13.927(018) 

13.591(016) 

1.335(032) 

14.395(067) 

1.874(127) 

0.00547 

0.025 

(19) 

2002el 


4 

16.153(025) 

16.179(020) 

1.375(025) 

17.381(044) 

1.961(139) 

0.02238 

0.085 

(20) 

2002er 


4 

14.266(012) 

14.108(012) 

1.253(010) 

15.070(037) 

1.984(079) 

0.00855 

0.157 

(21) 


R£FERENUE5^_^^fl^jLira et ~ iHajnuv et alj_|l^g^b^j_^f3j JAEavilla et_^J (4) l^imtzefl 

(5) lRj£h£nona etay^^ n^OsT T^^T 

(11) ; (14) 

( Isilyink^^t^Ll i200^ : (17) 

Krisciunas_et_^ 1 (19l lBenett^^!' m2004l) T(20‘) Wang, L. & Li, W.D. 2005, in preparation; (^ij^^^^^^e^ay 
i2004II . 

^The uncertainty in this quantity has been adjusted as discussed in Appendix|^ 
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0.6 0.8 1 1.2 1.4 1.6 

Aw 15 

Fig. 2.— Maximum-luminosity color E = vs. de¬ 

cline rate Amis calculated from our light curve fits. Closed squares 
(colored violet in electronic edition) show SNe from elliptical (E) 
and E/SO host galaxies; open squares (colored blue in electronic 
edition) show SNe from SO hosts. SNe from remaining (late-type) 
hosts are divided into those with normal color, shown by closed 
triangles (colored green in electronic edition); and those that are 
reddened, shown by closed circles (colored red in electronic edition). 
Open circles depict SNe that do not participate in the global fit. 


dust in the line of sight to the SN; and the photome¬ 
try references that supplied the data for our fits. All fit 
quantities have been AT-corrected. The tabulated uncer¬ 
tainties were propagated from the published photometry 
errors; uncertainties that are flagged in Tablets footnote 
were adjusted as discussed in Appendix IbI The origins 
of the SNe i n Table El are diverse , but they are draw n 
mai nly fromiHamuv et alJ lll996bll : iRiess et ’HI lll999ll : 
and iJhal (1200211 . The SN fliixes that were input to our LC 
and CMAGIC fits were taken from the compilation used 
by Conl05. 

All Imaxj magnitudes in Table [3 have 

b een corrected fo r MW d ust extinction using the dust map 
of iSchlegel et a, 1.1 (|1 99SI1 and the law of Card89. The cor¬ 
rections to and W™ used extinction coefficients 

Rb = 4.15 and Ry = 3.14, respectively, while, following 
equation 0 , the correction to used a coefficient 

that was smaller than Rb by the value of (3bv measured 
for the particular SN being corrected. More generally, 
quantities appearing anywhere in this paper that require 
correction for the effects of MW dust, whether they are 
numbers or symbols in equations, should be assumed al¬ 
ready to have been corrected for these effects. If neces¬ 
sary, our numerical corrections can be undone by using 
the values of (B — V) found in Table |3 

To help visualize the LC fit results. Fig. [3 displays a 
scatter plot of maximum-luminosity color E vs. Atoi 5 • 
As detailed in the caption, different symbols divide the 
sample by host-galaxy morphology. In 0the character¬ 
istics of the portion of Fig. |3in which SNe are clustered 
ut \E\ ^0-1 are explored. 


4. COMPARISONS OF Bmax TO Bbv 

Bmax can be compared to the CMAGIC output Bbv 
with substantially higher precision than is achieved by 
subjecting either of these luminosity measures to its own 
Hubble fit. Primarily this is true because Rmax — Bbv 
is immune to fluctuations in SN peculiar velocity that, 
in our range z < 0.1, are a major source of uncertainty 
in Hubble fits. (As well, Hmax — Bbv is blind to cos¬ 
mological variations in Hubble’s law.) After correction 
for observed SN color has been performed (as is im¬ 
plied thoughout this paper if the notation is absent), 
BmaK—BBv has the additional advantage that it depends 
only multiplicatively upon the effective color-correction 
coefficient 77., a phenomenological parameter that must 
separately be fit. As eqs. © and show, one may di¬ 
vide Hmax — Bbv by 77 to obtain the difference AA = 
£ — E between the CMAGIC color £ = {B^'^^ — B'^^y)/P bv 
and the maximum-luminosity color E = It 

is to the study of AA, which is completely independent 
of 77 as well as SN peculiar velocity, that this section is 
devoted. A drawback of studying /S.E is that any sys¬ 
tematic variation it may exhibit cannot, without more 
information, be apportioned uniquely between Hmax and 
Bbv- 

Because the observed SN color is influenced by intrin¬ 
sic SN color as well as intervening dust, and because 
Hmax and Bbv probe epochs at which intrinsic SN color 
might possibly be manifested differently, in principle AA 
might vary systematically with observed color itself. Af¬ 
ter correcting AA for systematic dependence on SN de¬ 
cline rate, to which we shall return shortly, we tested this 
possibility by fitting AE to a linear function of E. The 
additional free parameter allowed to diminish by only 
0.1, lending no support to this possibility. Accordingly, 
in the rest of this paper, we make no systematic color 
correction to AE itself, and, in results that are based on 
eqs. and ®, we draw no distinction between values 
of 77 used for calculating corrected Hmax as opposed to 
Bbv- 

4.1. Decline-Rate Dependence of Color Difference 

Figure 13 shows the dependence upon the decline-rate 
parameter Arn-is of the color difference AE. Though we 
were aware of the gentle nonlinearity in Amis of Hubble 
residuals for B and V already reported by Phil99, we 
were unprepared for the striking appearance of Fig. |3 
which has acquired the nickname “dogleg plot”. The 
data there are fit to a bilinear parametrization 

0 = AA — aA£;(Ami 5 — 1.1) — Amis ~ 1-1| -I- const , 

( 8 ) 

with a “kink” at the nominal centroid Amis = 1-1 mag. 
(Evidently, for fast decliners (Amis > 1-1) the decline- 
rate coefficient is uae + and for slow decliners it 

is Of AS — If only photometric errors on AE are 

taken into account, is unacceptably large (99.5 for 
53 degrees of freedom (dof)). Adding in quadrature to 
the photometric errors a modest rms noise {i.e. intrin¬ 
sic scatter) term n^E = 0.029 mag reduces to its 
most probable value of 51, where its probability density 
is larger by the factor exp (7.2). The same likelihood ra¬ 
tio would be associated with a gaussian fluctuation of 
\/14.4 = 3.8 standard deviations (3.8cr). Therefore, in 
the difference between CMAGIC and maximum-luminosity 
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-0.4 -0.2 0 0.2 0.4 

A£ uncorrected for Am 15 dependence 

Fig. 3.— Decline rate Amis (m^g) the difference AE (mag) 
between CMAGIC and maximum-luminosity colors. In this plot the 
abscissa is uncorrected for decline-rate dependence. The fit line 
is the function in equation e. SNe marked by open circles are 
omitted from the fit. 


B — V color, intrinsic noise clearly is present at a modest 
level (0.029 ±0.007 mag). 

Only after noise is added to bring down to a prob¬ 
able level (x^/dof —> 1 as dof ^ 00 ) can the errors on 
the fit parameters be trusted. These parameters and 
their errors are shown in the upper portion of Table 13 
Various alternate conditions for the fit are exhibited in 
the lower portion of the same table, together with the 
corresponding Ax^ penalty that is assessed when each 
condition is enforced. The AE data definitely reject the 
hypothesis of no Amis dependence at all (Ay^ = 86.1 
for 2 dof). Fitting the data to the bilinear parametriza- 
tion from equation both a nonzero linear coefficient 
o^AE (7.5a) and a nonzero kink coefficient (5.8a) are 
required. Therefore, in AE the presence of a nonlinear 
decline-rate dependence is firmly established. 

Especially because the bilinear parametrization is non¬ 
standard, it is natural to explore alternatives. Obvi¬ 
ously, the first derivative with respect to Amis of any 
parametrization ideally should be continuous. For the 
bilinear case this could be accomplished by substituting 
a small curved segment for the corner at Amis = 1-1. 
Here we omit this feature because the data near the cor¬ 
ner are too sparse to determine the curvature there. The 
most common nonlinear form is a quadratic. As recorded 
in the final row of Table|3 a Ay^ penalty of ±7.5 is paid 
for substituting a quadratic function of Amis for the bi¬ 
linear form. Statistically, with respect to the bilinear 


TABLE 3 

Results of bilinear fit to data in 
Fig. El 


Result of AE fit 

Fit value 

riAE 

0.029 ± 0.007 

aAE 

0.23 ±0.03 

“ab 

-0.30 ±0.05 

X^/ dof 

51.0 / 53 

Fit conditions^ 

Cl 

<1 

default (as above) 

0 

riAE = 0 

±48.5 

“AB = 0 

±57.1 

“ab = d 

±33.6 

(^AE — ^ae ~ ^ 
quadratic (not bilinear) 

±86.1 

±7.5 


^ Other parameters are allowed to vary. 


form, this sets the betting odds against the quadratic 
at 43:1. Another factor weighing against a quadratic 
parametrization is its asymptotic behavior: the data in 
Fig.|3do not appear to prefer a steeper slope in Amis as 
the extremes in that parameter are approached. For all 
further analysis in this paper we continue to parametrize 
the Amis dependence by a bilinear function like that in 
equation ®. 

With intrinsic noise and nonlinear decline-rate depen¬ 
dence established as basic features of — Ssy by the 
results in Fig. |3 and Table |3 the key issue not yet ad¬ 
dressed is raised: how is responsibility for these features 
apportioned between Bmax and Bbv7 This statistically 
challenging issue must be attacked using Hubble fits, to 
which we now turn. 

5. HUBBLE FIT TO AVERAGE OF Hmax AND Bbv 

With a relatively precise fit to the difference of B^aax 
and Bbv now in hand, we choose as the input to our 
primary Hubble fit the orthogonal linear combination of 
these two magnitudes, i.e. t heir u nweighted average (B). 
As will be demonstrated in ^5.31 the sum of y^ for both 
of these fits is equivalent to the y^ for a global fit to all of 
the Hmax and Bbv data. Although this is the simplest 
and most straightforward type of global fit, it is not the 
only possible approach: for example, one could set aside 
the fit to Bmax — Bbv and substitute simultaneous Hub¬ 
ble hts to H„iax and Bbv- That approach would have re¬ 
quired taking careful account of the {Hmax, Bbv} covari¬ 
ance matrix for each SN, including the large off-diagonal 
term due to peculiar velocity. Using that approach, the 
striking decline-rate dependence shown in Fig. 13 could 
have been reproduced, but the data supporting it would 
have been buried in the details of those matrices. 

All Hubble hts performed in this paper assume a hat 
universe with the ratio Dm of mass to critical densities 
set at 0.3. Over the range in redshift of the SNe we 
consider, z < 0.1, the results of these hts are insensi¬ 
tive to variations in Dm and Oa that reasonably would 
be allowed by current world data. Also, for every Hub¬ 
ble ht that is reported here, as well as for every ht to 
AE = 0, a constant overall magnitude offset is included 
that always is allowed to assume its own best ht value, 
independent of the offset determined by any other ht. 
(For Hubble hts, this offset is essentially the absolute SN 
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magnitude combined with the Hubble constant.) The 
best fit values of these offsets are left unreported in this 
paper since they are irrelevant to its goals. Because only 
linear color corrections are used here, any constant offset 
in color plays the same role as a constant overall offset 
and therefore leaves unaffected any fit results that we 
report. Accordingly this paper does not distinguish be¬ 
tween a color B — V and its excess with respect to an 
expected color. 

Apart from the assumptions implicit in using particu¬ 
lar B and V light-curve templates, and in exploiting the 
observed linear relationship between B magnitude and 
B — V color in the CMAGIC region, all analysis in this pa¬ 
per makes no prior assumption. Specifically, because the 
effects of extinction by host-galactic dust and of intrinsic 
SN color variation are not yet well characterized, we ac¬ 
cept and use the observed B — V color without applying 
to it any correction based e.g. on a Bayesian prior. The 
statistical error analysis applied here is frequentist: in 
a multiparameter fit, when a single parameter is varied 
and the change in likelihood is recorded, all other pa¬ 
rameters (including the offset) are allowed also to vary 
for best compatibility with that excursion. (If instead we 
had chosen to integrate over a likelihood surface, implic¬ 
itly we would have had to apply a Bayesian prior.) Thus 
our calculated errors are never assumed to be symmetric; 
when an average error is reported, the calculated errors 
lie within ±5% of it. 

5.1. Correcting (B) for Color and Decline Rate 

Because AB = Hmax — Bbv is proportional to AE 
(eqs. (|21) and ®), the bilinear dependence of AE upon 
Atoi 5 should map into a bilinear Ami 5 -dependent cor¬ 
rection to Bmax, Bbv, or both. The Hubble fit described 
here allows for both. Allowing for corrections for both 
color and decline-rate dependence, the more general form 
of eqs. © and © becomes 

Hmax = Blf;Zr -TIE- I?(Ami 5 ; ^max, Omax) (9) 

Bbv = BI^Zl - TZS - 'D{Ami5-,aBv, oc'bv) , ( 10 ) 
where the function D{Ami^]a,a') is shorthand for the 
bilinear decline-rate dependence introduced in equation 

T>(Ami 5 ; a, ct') = a(ATOi 5 -l.l)-|-a'|Ami 5 -l.l| . (11) 

Apart from the “kink” term in the corrections 

described by eq u ation m are equivalent to those intro¬ 
duced bv iTrionl lll998ll . The coefficients aBV and a'^y 
for Bbv are not the same as the coefficients amax and 
Q^max ^max; rather, the two pairs of coefficients are 
linked by the decline-rate corrections to AE = £ — E 
given by equation ®. Straightforwardly, the average 
(B) of Hmax and Bbv is corrected for color and decline 
rate by a similar relation: 

(B) = B]^:^ - n{E) - I?(Ami5; ^(b), ^(b)) , (12) 

where {E) = (A -|- £)/2 and the decline-rate coefficients 
oti^B) and are equal to the average of those for 

Hmax and Bbv- Table 0 provides a full list of the bilin¬ 
ear decline-rate coefficients and their interrelationships. 
There, in terms of the color coefficient TZ and of the 
decline-rate coefficients for correcting AE and Bmax, the 
last column supplies the decline-rate coefficients used to 
correct AH, (H), and Bbv- 


TABLE 4 

Definitions of bilinear decline-rate correction 

COEFFICIENTS USED IN THE FIT TO AE = 0; USED FOR 
correcting AB; AND USED IN HUBBLE FITS TO 

BmaK — n, — n, and Bbv — n- 


Coef¬ 

ficient 

Multiplying 

To 

correct 

Is equivalent 
to 

OtAE 

(Amis 

-1.1) 

AE 


“as 

|Ami5 

-1.1| 

AE 


“AS 

(Amis 

-1.1) 

AB 

7?.«ab 

“as 

|Amis 

-1.1| 

AB 

^“ab 

Ciimax 

(Amis 

-1.1) 

-^max 

^max 

|Amis 

-1.1| 

-^max 

71 

Q^max 

*^max ~2^AE 

“(B) 

(Amis 

-1.1) 

(B) 

“(B) 

|Ami6 

-1.1| 

(B) 

“SV 

(Amis 

-1.1) 

Bbv 

<^max — '^O'AE 
*^max '^^AE 

“sv 

|Amis 

-1.1| 

Bbv 


For completeness, this section concludes with a brief 
discussion of the “effective color” that multiplies the 
color-correction coefficient TZ in eqs. itTIIIl and (tT^ . This 
detail will be relevant only for visualizing the quality of 
the fit that determines TZ. In eqs. itTIIIl . and itT^ . the 
coefficients TZ are all the same, but they multiply colors 
- E, £, and (E) - that obviously are different. Less 
obvious is the fact that the right-hand term in each of 
eqs. ns and m itself is a linear function of TZ, through 

the dependence of a^BV '^(b) ^ that is shown 

in the last four rows of Table 0 Using the information 
tabulated there to reexpress cx^bv "^(b) terms of 

ccmax and and collecting into H®® the sum of the 

coefficients of TZ in equation m and in equation 
one obtains the effective colors 

Bbv- = £ — D{Ami^', a/^E, (x'^b) (13) 

(B) ■- = (E) - X>(Ami5; ^oab, I^ab) ■ (14) 

When it is uncorrected for color, the Hubble residual 
e.g. of (B) should exhibit a linear dependence on H®®^ 
with slope TZ. 

5.2. Propagating Intrinsic Noise into (B) 

Leaving aside sources of noise that are common to both 
Hniax and Bbv, and introducing a mild assumption that 
is discussed fully in Appendix IeI the rms intrinsic noise 
Bab that is measured in the comparison of Hmax to Bbv 
propagates into (H) as n(^B) = nABl‘^- Therefore, using 
the factor TZ to relate the intrinsic noise in AB to that in 
AE, the added rms intrinsic noise bab = 0.029 ± 0.007 
mag that was required by the fit to AE = 0 (Table |2l 
propagates into an added rms intrinsic noise term ni^B) = 
77-bab/ 2 in the average (H) of Hmax and Bbv- 

As is also discussed in Appendix|El intrinsic noise that 
is common to both Hmax and Bbv may be present, to¬ 
gether with common-mode noise from the effects of pecu¬ 
liar velocity along the line of sight. In principle, measur¬ 
ing the redshift dependence of the common-mode noise 
would allow the effects of peculiar velocity to be resolved 
from those of other sources. However, the SN sample 
statistics available to us do not permit a clean separa¬ 
tion. Instead we resort to parametrizing the sum of all 
sources of common-mode noise as the manifestation of a 
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Pull (A£) 

Fig. 4.— Pull (= residual/error) from Hubble fit to (B) vs. pull 
from fit to AE = 0. SNe that lie within the circle (filled points) 
contribute Ax^ < 8 to the global and are included in the fits. 

single effective rms SN velocity tieg relative to the Hub¬ 
ble flow. In 45.01 we report the value of VeS determined 
from the data, from which we deduce an upper limit on 
the true rms peculiar velocity. 

5.3. Equivalent Global Fit and Fitted Sample 

Figure 0 shows for each SN the pull (defined as the 
residual divided by the error on that residual) for the 
Hubble fit to (B) vs. the pull for the fit to AE = 0 
described in SH The insignificant correlation between 
the two pulls, measured by the Pearson correlation co¬ 
efficient |R| < 0.05, confirms that, as expected, the two 
fits are orthogonal. Therefore the sum of for both 
fits is equivalent to the for a global fit to the param¬ 
eters aAE, c('ae^ ^ (plus offsets). Note 

that the decline rate parameters for Hmax and Bbv are 
also fully specified by the global fit (Table 0. One SN 
contributes two degrees of freedom to the global ht: one 
dof for the Hubble fit to (B) and one dof for the fit to 
AE = 0. Its contribution to the global y^ is equal to the 
sum of the squares of the same pulls exhibited for that 
SN in Fig. ^ 

If it lies outside the circle of radius \/8 displayed in 
Fig.EI one SN contributes to the global y^ an increment 
Ay^ > 8 for 2 dof. If the pulls were distributed according 
to gaussian statistics, for our sample size an average of 
one SN would be expected to lie outside this circle. In 
fact five SNe lie outside: four mainly because of a large 
I pull I in the Hubble fit, and one primarily because of a 
large |pull| in the fit to AE = 0. As previously noted in 
m we exclude from fits all five of these SNe, leaving a 
fitted sample of 56 SNe. Using open-circle symbols we 
continue to exhibit these outliers in plots. 

5.4. Decline-Rate Correction to Hubble Fit for {B) 

In Fig.Elthe decline rate Atoi 5 is plotted vs. the Hub¬ 
ble residual for {B). In this plot (B) is fully corrected 



({B}- jii) uncorrected for Am j 5 dependence 

Fig. 5. — Decline rate Amis Hubble residual ({B) — /j.), 
which here is uncorrected for its dependence on Amis. The line 
is a bilinear function of Amis from the global fit, from which SNe 
marked by open circles are omitted. 


for the effects of color, as will be discussed in 45.51 but 
it is left uncorrected for decline rate. The line is the 
correction determined by the global fit and described by 
equation US- The best-fit kink at Amis = 1-1 is less 
pronounced than that observed for AE (Fig. |3). Be¬ 
cause the decline-rate dependence is much smoother for 
the sum of Umax and Bbv than for their difference, the 
kink in decline-rate dependence should (and does) ap¬ 
pear with opposite sign in the Hubble residuals f or Hm ax 
and Bbv', this point will be discussed further in ED In 
Table El the top two rows list the globally fit coefficients 
of Amis and |Amis that correct the Hubble residuals 
of {B) via equation 1I2II ; these parameters determine the 
line in Fig. El 

5.5. Color Correction to Hubble Fit for (B) 

Figure El exhibits the color correction to {B). On the 
ordinate is plotted the effective color from equation 
1I4II : this is the coefficient of TZ in equation 1I2II . As in 
Fig. El the horizontal coordinate is the Hubble residual 
for {B ); however, in this plot (B) is fully corrected for the 
effects of decline rate, but it is left uncorrected for color. 
Again, the line is determined by the global fit; in this 
case the correction is purely linear, with best fit slope 
TZ = 2.59 ± 0.24. This result is inconsistent either with 
a null correction or with a coefficient typical of Milky 
Way dust: for the alternate values 7^ = 0 or 7^ = 4.1, 
respectively, y^ increases by >115 or >32. 

When the color is measured with less than perfect 
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TABLE 5 

Bilinear decline-rate correction 
COEFFICIENTS EOR (B) , B max} AND 
BbV from the global FIT. 


Result of global fit 

Fit value 

«(S) 

0.48 ± 0.09 

rv' 

(-B) 

-0.14 ± 0.18 

^max 

0.78 ± 0.09 

Q/^ 

*-*max 

-0.54 ± 0.19 

aBV 

0.18 ± 0.12 

^BV 

0.25 ± 0.21 

»max — CiBV 

0.60 ± 0.09 

Ol' - Ol' 

^max ^BV 

-0.79 ± 0.16 

Fit conditions^ 

Ax^ 

default (as above) 

0 

“max = 0 

-1-7.2 

ttmax + D^max “ ^ 

-1-2.3 

“sv = “sv = 0 

-1-9.2 

Qmax ^BV — 0^ 

+57.1 

^max ~ ^'bV ~ 

+33.6 


® Other parameters are allowed to vary. 
'^From Table 13 



{{B)- ju) uncorrected for color dependence 

Fig. 6. — Effective color (= coefficient of TZ, defined in 

equation 1141 1 vs. residual from Hubble fit to {B) when the latter 
is left uncorrected by terms that are proportional to TZ. The line is 
a linear function of the ordinate determined by the global fit, from 
which SNe marked by open circles are omitted. 


precision, the slope TZ will be underestimated. In this 
analysis, the effects of finite color resolution were un¬ 
folded by two different methods described in Appendix 
IfI These methods yielded corrections to TZ of -1-3.8% 
and -1-7.0%, respectively. For comparison to other de¬ 
terminations of TZ, which are likely to involve color res¬ 
olutions similar to ours, we prefer to quote the uncor¬ 
rected value TZ — 2.59 ± 0.24. On the other hand, if our 
value of TZ is to be compared to other determinations for 
which the effects of color resolution already have been 
unfolded, we estimate our own correction for those ef¬ 
fects to be given by the average of our two methods: 
/S.TZ/TZ — -1-0.054 ± 0.027, where a 50% systematic error 
is assigned to the correction. 

As an essential check, we have used the results of pseu¬ 
doexperiments based on two independent Monte Carlo 
simulations to study the distributions both of pulls and 
residuals in TZ and other fit parameters. In these pseudo¬ 
experiments, minimizing global to obtain TT-et does re¬ 
produce on average the color-correction coefficient 72.true 
that was input to the simulation; measuring the variation 
of with TZf^t does predict the dispersion of TZf^t — T^true 
that is measured. As is conventional, our global fit min¬ 
imizes = 5^ Ti 5, where i5 is a vector of differences and 
Ti is the information (inverse covariance) matrix. This 
minimization problem belongs to the class in which both 
5 and Ti depend on the fit parameters; the latter depen¬ 
dence occurs because the uncertainties in Hubble resid¬ 
uals grow e.g. as TZ increases. For these Monte Carlo 
checks to be satisfied, we found it essential to allow Ti. 
to vary continuously as the fit parameters were being 
optimized. 

We also tested an iterative stepwise fitting technique in 
which was minimized with Ti held fixed, then Ti was 
recalculated using the newly updated parameters, etc., 
until the iteration converged. This caused our Monte 
Carlo checks to fail badly. If applied to the real SN sam¬ 
ple analyzed here, the stepwise technique would yield a 
value of TZ that is 10.2% below that obtained with the val¬ 
idated fitting procedure. Also, the measures of Hubble- 
line deviation to be discussed in H6.31 would be marginally 
reduced (by «3%). This is understandable: slightly un¬ 
dercorrecting e.g. for color variation can marginally im¬ 
prove the Hubble dispersion if the uncertainty associated 
with that correction can be reduced. Nevertheless, all 
results presented in this paper are based on the Monte 
Carlo-validated fitting technique described in the previ¬ 
ous paragraph, even at the expense of marginally larger 
Hubble dispersions. We made this choice not only be¬ 
cause the color- and decline-rate correction coefficients 
are physical properties of the SNe that we wish to re¬ 
port accurately, but also because in future campaigns, as 
the photometric errors e.g. in color become relatively less 
important, the correction coefficients which minimize the 
Hubble dispersion will approach their true values. 

5.6. Upper Limit on Peculiar Velocity 

As described in a single effective rms SN velocity 
Ueff is used to parametrize the effects both of peculiar 
velocity Upec along the line of sight relative to the Hubble 
flow, and of any other non-photometric noise component 
that is common to Hmax and Bbv- In the Hubble fit for 
{B), Veff is adjusted to maximize the probability density 
of x^, yielding the best fit value UeS = 382^5° km s“^. As 
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( 5 jnax“ /O uncorrected for Am 15 dependence 

Fig. 7.— Decline rate Amis Hubble residual (Smax — m), 
which here is uncorrected for its dependence on Amis. The line is 
the bilinear function of Amis from the global fit, from which SNe 
marked by open circles are omitted. 

noted in measuring UeS yields an upper bound on 
SN peculiar velocity. The corresponding 95% confidence 
upper limit is Upec < 486 km s“^. 

In the literature, estimates are available for the pe¬ 
culiar velocities of individual SNe in the Local Group; 
applying corrections based on those estimates would im¬ 
prove significantly our Hubble fits. However, our aim 
here is to use local SNe to characterize the peculiar ve¬ 
locities of distant SNe, for which such estimates generally 
are unavailable. Therefore the results quoted above for 
Uefi and Upec intentionally were left uncorrected for the 
estimated velocities, relative to the Hubble flow, of local 
SNe. 

6. RESULTS FOR Hmax AND Bbv 

As emphasized in >15.hi the fit to AE = 0 and the Hub¬ 
ble fit to {B) together comprise a global fit that fully 
determines the decline-rate and color corrections to be 
applied to Bmax and Bbv as well as to (B) and AE. In 
principle, no separate Hubble fits to Hmax or Bbv are re¬ 
quired. Figure[7|plots the decline rate Amis vs. the Hub¬ 
ble residual for Hmax- The corresponding decline-rate 
coefficients for Hmax (and also for Bbv) are supplied in 
Table El There the symbols a and a' refer, respectively, 
to coefficients of (Amis — 1-1) and |Amis — 1-1|; they are 
fully defined in Tabled] For alternate fits in which one 
or more decline-rate coefficients or combinations thereof 
are forced to vanish, the lower part of Table El exhibits 
the corresponding Ax^ penalty. 


Strovink et al. 

As a check, independent Hubble fits to each of i?max 
and Bbv were performed, allowing TZ and both bilinear 
decline-rate coefficients to seek optimum values indepen¬ 
dent of the global fit. Relative to values of obtained 
for Hubble hts to Hmax and Bbv using parameters deter¬ 
mined by the global ht, values of for these independent 
fits were smaller by less than 0.8 and 2.5, respectively. 
This confirms that the global fit is a reasonable descrip¬ 
tion of both Hmax and Bbv- 


6.1. Decline-Rate Corrections to Hmax and Bbv 

As a guide to the variety of numerical results from 
the global ht that are included in Table El first we high¬ 
light those features which are firmly established by the 
data. In particular, the results exhibited in rows 7 and 
8 of Table El derived from results of the fit to AE = 0 
summarized in Table El emphasize that Hmax and Bbv 
cannot share the same decline-rate coefficients. Because 
a kink in decline-rate correction at Amis = 1-1 is a fea¬ 
ture of Type la SNe that may not be familiar to every 
reader, here we focus particularly on the coefficients a' of 
I Amis — 1.1 1 . Because the tabulated values of aj^ax and 
a'gy are nonzero, respectively, only by 2.7 and 1.2 stan¬ 
dard deviations, a reader might imagine that the data 
marginally would allow both to vanish. However, be¬ 
cause the difference between those coefficients is nonzero 
at the 5.8tT level, at least one of Hmax and Bbv must be 
described by a nonzero kink. Other firmly established 
features include the need for a conventional decline-rate- 
correction term that is linear in Amis for both (B) and 
7?max (rows 1 and 3 of the table). 

Conversely, the data may allow certain coefficients in 
Table El or combinations thereof, to vanish. As one ex¬ 
ample, the result for o.'(^b) found in the second row of 
the table allows the decline-rate correction to (B) to be 
purely linear; that would correspond to the case in which 
the fit line in Fig. El is straight. As a second example, 
fast-declining SNe (with Amis > 1-1) are corrected by 
a coefficient of (Amis — 1-1) that is equal to a -I- a' . 
As shown for corrections to Hmax iu the fourth-from-last 
row of Table El c«max + Q^max is only 1.5 standard devia¬ 
tions from zero. Therefore the data do not conclusively 
reject the proposition that fast-declining SNe require no 
decline-rate correction whatsoever to Hmax; that would 
correspond to the case in which the top part of the “dog¬ 
leg” line in Fig. [3 is exactly vertical. 

Regarding the CMAGIC luminosity Bbvi the third- 
from-last row of Table El examines the proposition that 
Bbv requires no decline-rate correction for any SN: when 
two dof are removed by nulling both ubv and 

grows by 9.2 (confidence level 0.01). Evidently the 
data disfavor this proposition. More restrictively, slow- 
declining SNe (with Amis < 1-1) are corrected by a 
(Amis —1.1) coefficient equal to a—a'] the results in rows 
5 and 6 of Table El clearly permit ubv — ch'bv to vanish. 
Therefore the data do allow slow-declining SNe to remain 
uncorrected for decline rate if analyzed using the CMAGIC 
method. For color as well, some of these SNe might ac¬ 
ceptably remain uncorrected: “bluish” SNe (with B — V 
color £ < 0.1) have a color correction to Bbv that is 
no greater than {TZ — Pbv)^ ~ (2.6 — 2) (0.1) « 0.06 
magnitudes. 
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TABLE 6 

Consequences oe allocating the intrinsic noise in AB 

ALL TO BbV, all TO B max? OR. EQUALLY TO BOTH. 


Quantity Noise allocation 

all to Bbv equal all to Bmax 


Conditions and Hubble fit results: 


Noise in Bmax (mag) 

0 

0.052 

0.074'* 

Noise in Bqy (mag) 

0.074'* 

0.052 

0 

(Bmax residuals) 

43.2 

39.3 

36.5 

{Bbv residuals) 

57.8 

66.1 

80.2 

Relative likelihood 

1 

1/3.3 

1/44 

^For {Bmax — Bbv) the fit noise is ' 

0.074 ± 0.019 

mag. 


6.2. Allocating Intrinsic Noise to Rmax or Bbv 

As explained in ^5.21 the intrinsic noise term n^E that 
is required by the fit to NE = 0 propagates into a noise 
term m^s) = in the Hubble fit for {B). Table El 

presents the for the Hubble fit to Hmax, the for the 
Hubble fit to Bbv^ and the relative likelihood (propor¬ 
tional to the product of the two densities) for three 
scenarios: (a) all of ni^B) is allocated to Bbv] (b) l/-\/2 

of ni^B) is allocated to Bbv and 1/a/2 to Hmaxi (c) all of 
ni^B) is allocated to Hmax- All of the fits that contribute 
to the table use the parameters of the global fit. Sce¬ 
nario (a) is favored, so it is the basis for the global fit 
results that are presented in this paper. Scenarios (b) 
and (c) are progressively less likely, by factors of 3.3 and 
44, respectively; evidently scenario (b) is not strongly 
disfavored relative to (a). Although the data do not con¬ 
clusively reject scenario (c), it is unlikely that all of the 
noise in the difference between Hmax and Bbv is due to 
noise that is confined to Bmax- 

6.3. Quantifying Deviations about the Hubble Line 

From the distribution of a property m of any statisti¬ 
cally limited sample, estimating the variance of the mea¬ 
sured variance in m of that sample requires no assump¬ 
tion about the shape of the true m distribution. In part 
for that reason, the unweighted rms Hubble residual a 
(the square root of the variance) is a useful standard 
measure of SN uniformity. A drawback of the rms is 
that, if included in the sample, SNe with unusually large 
uncertainties in m can dominate this measure, pushing 
the rms upward. This shortcoming may be addressed by 
use of the weighted rms (wrms), denoted by a^- The 
square of is the mean square deviation of m from the 
function to which it is fit, weighted by the inverse square 
uncertainty assigned to m for each SN. (Equations defin¬ 
ing both cr and are provided in Appendix 1X1 1 These 
definitions take proper account of the reduction in dof 
due to the number of variable parameters in the fit func¬ 
tion. 

While it is invariant to the overall scale of the uncer¬ 
tainties assigned to m, the wrms is fragile if the individ¬ 
ual uncertainties are assigned inaccurately. As a patho¬ 
logical example, if 100 identical SNe are observed with 
identical methods, and if the uncertainty in m of the first 
SN erroneously is set to 1% of the value assigned to the 
rest, the wrms underestimates the true rms typically by 
a factor of >7. 


TABLE 7 

Weighted rms deviation, rms deviation, and core 
DEVIATION (mag) EROM THE HUBBLE LINE OVER THE 
RANGE 0.015 < Z < 0.1 FOR {B), Bmax, AND Bgy USING 
THE GLOBAL EIT AND EOR Bgv USING A Ssi/-ONLY EIT. 


Hubble 

deviation 

of 

Type 
of fit 

Weight¬ 
ed rms 

rms 

1.483 X 
median 

1 residual! 

(B> 

global 

0.129 

0.152±0.018 

0.116 

-^max 

global 

0.130 

0.152±0.015 

0.144 

Bbv 

global 

0.155 

0.176±0.018 

0.158 

Bbv 

Bbv only 

0.132 

0.144±0.015 

0.120 


Both the rms and the wrms are highly sensitive to the 
presence or absence of a few SNe far out on the tails of the 
distribution. This sensitivity is addressed here by intro¬ 
ducing as a third measure of SN uniformity the median 
absolute residual meddm^ — /i(zi)|), or “core deviation”. 
We scale this quantity by the factor 1.483 so that it is 
equivalent to the rms for a pure gaussian distribution in 
the limit N/Ndo{ 1- 

6.4. Hubble-line Deviations of (B) , Bmax, and Bbv- 

The first three rows of Table [7| list the wrms, rms, 
and core deviations from the Hubble line for (B), Bmax, 
and Bbv^ using the single set of correction coefficients 
obtained from the global fit. To suppress the effects of 
peculiar velocity, only SNe with z > 0.015 are included 
in these averages (Fig. El exhibits the sample’s z distri¬ 
bution). We estimate the uncertainties in values of the 
weighted rms and core deviation to be similar to the rel¬ 
atively well-defined uncertainties in the values of the rms 
that appear in the table. That these uncertainties are as 
large as their estimates is emphasized by the final row of 
Tabled which shows the result of substituting for the 
output of the global fit the output of a Hubble fit that 
was made only to Bbv- Residuals from the Bsy-only fit 
are shown in Fig.E) In this fit, after adding 0.053 mag of 
noise in quadrature to the photometric errors in Bbv i we 
rejected as outliers five SNe chosen to minimize for the 
Bsy-only fit rather than for the global fit. The result¬ 
ing probability density was maximal. Two SNe that 
had been excluded from the global fit, SN1992bp and 
SN1994ae (Fig.®, were replaced as outliers by SN1995ac 
and SN2001V. As is apparent, substituting the B^v^-only 
for the global fit perturbed the Hubble-line deviations for 
Bbv by amounts comparable to other differences found 
in the table. Weighing all of the information in Table 
d it is evident that, with the available statistics, the 
Hubble-line deviations for all three measures of SN lumi¬ 
nosity - (B), Bmax, and Bbv ~ are comparable and of 
order «0.14 mag, and are not easily distinguished from 
each other. Note that these deviations arise from the ef¬ 
fects of photometric error and peculiar velocity as well 
as intrinsic noise. 

7. CHARACTERISTICS OF THE BLUE RIDGE 

SNe that are considered to be part of the “blue ridge” 
are grouped in a region roughly defined as |Bmax — 
KnaxI 0-1; the word “ridge” refers to the contour of 
this group on a plot of color vs. decline rate, as in Fig.d 
SNe in the blue ridge are thought l argely to be free of 
extinction due to host-galactic dust llHatano et alJll99'^ 
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{Bgy - ju) fully corrected 

Fig. 8. — Logio(redshift z) vs. OMAGIC Hubble residual [Bgy — 
fi). The points are residuals from the Esv-only fit. Five SNe 
excluded as outliers are suppressed. The curve shows the itlcr 
uncertainty in residual associated with the assumed rms line-of- 
sight peculiar velocity of 400 km/s relative to the Hubble flow. 
Accounting for the dof removed by the corrections, the unweighted 
rms deviation from the fit (straight line) of the points with z > 
0.015 is 0.144 ± 0.015 mag (Table|7J. 

ICommin3l2004D . If this is so, sufficiently well-observed 
variation of SN color within the ridge area could reveal 
the distribution of intrinsic SN color. Here we examine 
whether such observed color is correlated with decline 
rate Arn-is, and whether SN luminosities require correc¬ 
tion for its variation. 

7.1. Decline-Rate Dependence of Ridgeline 

We have searched for evidence of a linear dependence 
upon Atoi 5 of the blue ridgeline by making an unbinned 
maximum-likelihood fit of the data in Fig.|21to a “ridge” 
color-distribution function (whose most probable color 
varies as A{Ami 5 — 1.1), where yf is a fit constant), 
plus a Amia-independent “tail” color-distribution func¬ 
tion that is nonzero only on the red side of the ridge 
peak. (For the narrow purpose of this fit, the advan¬ 
tage of this simple color-distribu tion function o v er a m ore 
realistic model, such as that of lHatano et al.l lll998D . is 
that SNe far outside the ridgeline have smaller influence 
on the fit results.) Within statistics {cr{A) « 0.06), we 
find no evidence for a nonzero value of the coefficient A. 
Because our blue ridge is consistent with being Amia- 
independent, the only variations with decline rate of our 
corrections for SN color are those represented by eqs. d 

[TI| . 


7.2. Do SNe in the Blue Ridge Require Color 
Correction ? 

If SNe that are clustered in the blue ridge are largely 
unextinguished by host-galactic dust, so that they ex¬ 
hibit mainly a spread in intrinsic color, any necessary 
color correction to their luminosities could depart sub¬ 
stantially from corrections that normally are applied to 
extinction by dust. In particular, the data might require 
little or no correction to luminosity for intrinsic color 
variation. If this simple picture were correct, restricting 
analysis to SNe that are clustered in the blue ridge could 
substitute for implementing a color correction. 

To test this simple model, we need to identify those 
SNe which truly belong to the blue ridge group. To 
perform this identification while avoiding any chance of 
prejudice against the model, adventurously we allow the 
Hubble residuals of the SNe, as well as their colors, to 
play a role in the process of SN selection. For simplicity 
and for ready comparison with other work, in this section 
we consider only the Hubble residuals for i?max, and, in 
Hubble fits to Umax, we allow all extraneous parameters 
including decline-rate coefficients to assume their best fit 
values independent of the results of any other fit. 

Our scheme for identifying a group of SNe that have 
a relatively good chance of obeying the simple model 
described above is illustrated in the scatter plot shown 
in Fig.|^ There, for each SN, the abscissa is its pull in the 
Hubble fit to Umax assuming no color correction, i.e.TZ = 
0. The ordinate is a “color pull” whose numerator is the 
excess of Bmax — Vmax above a typical value of 0.02 mag, 
divided by the sum in quadrature of the photometric 
error and a characteristic intrinsic color rms of 0.03 mag 
taken from Phil99. As in Fig. 01 SNe that lie within 
the circle have a sum of pulls squared that is less than 
8. These 41 SNe comprise the sample that we select 
to have a relatively good chance of obeying the “TZ = 
0” simple model. Among the SNe rejected by this cut 
are SN19930, which is hosted by an elliptical galaxy; 
and SN1992bh and SN1995al, which belong to a range of 
observed color that we do find to be characteristic of the 
observed colors of SNe that are hosted by E, E/SO, and 
SO galaxies. 

It may appear naive to accept as many as 41 of the 
56-SN fitted sample (Table ^1 into a sample to be tested 
for freedom from extinction. However, the mean and rms 
Bmax — Knax for all 41 of these SNe (0.022 and 0.061, re¬ 
spectively) are comparable to those (0.022 and 0.062) of 
the sample of 13 SNe hosted by E, E/SO, or SO galaxies, 
whose reddening would be expected to be less significant 
than for the balance of SNe that occurred in spiral galax¬ 
ies. 

The final step in testing the simple “TZ = 0” model is 
to subject this 41-SN subsample to two alternative Hub¬ 
ble fits: one with TZ — 0 and the other with TZ = 2.59, 
the result of the global fit described in The details 
of the SN selection and fit results are shown in Table 
IHl In the top half of the table, the color and Hubble 
pulls are those described above. In the bottom half, the 
values of are for Hubble fits to with TZ fixed to 
either of the two indicated values. Even for the restricted 
subsample of SNe that is fit, for the TZ = 2.59 fit is 
substantially smaller than for 77 = 0, by a difference of 
13.8. Therefore, with 99.9% confidence, the data reject 
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Pull (color-uncorrected 5 njax - M ) 

Fig. 9.— “Color pull” (defined in the text) vs. pull in a color- 
uncorrected Hubble fit to iJmax, for SNe in five classes depicted by 
the same symbols (and colors in electronic edition) that are used 
in Fig. 1^ SNe that are clustered within the circle of radius \/8 
are identified as having a relatively good chance of not requiring a 
correction for color variation. 


TABLE 8 

Details of SN selection and results of fits that reject a 

SIMPLE MODEL IN WHICH SNe THAT ARE CLUSTERED WITHIN THE 
BLUE RIDGE REQUIRE NO GOLOR CORRECTION. 


Quantity 

Inside cluster 

Outside 

cluster 

PulF (color) -1- pulF (Hubble) 

< 8 

> 8 

SNe with E, E/SO, or SO hosts 

12 

1 

SNe with late-type hosts 

29 

14 

Conditions and results for Hubble fits to B 

max* 

IZ (fixed in fit to 41 SNe) 

2.59=^ 0 


(of fixed-7^ fit to 41 SNe) 

31.8 45.6 


Relative likelihood 

1 0.001 



®Prom global fit to 56 SNe, 7?. = 2.59 it 0.24. 


the hypothesis that SNe in the blue ridge group require 
no color correction. Whether this correction is for intrin¬ 
sic color or for dust, it is needed. Had we selected the 
test sample less aggressively - for example, had we de¬ 
clined to include among the selection criteria the sizes of 
the TZ = 0 Hubble residuals themselves - this conclusion 
would only be strengthened. 

7.3. Color Corrections to SNe with E, E/SO, or SO 

Hosts 


Of the 41-SN subsample chosen for the test just de¬ 
scribed, 12 are hosted by elliptical, E/SO, or SO galax¬ 
ies. A thirteenth, SN19930, has an elliptical host but 
was rejected from that subsample because, if uncorrected 
for color, it exhibited a large 7^ = 0 Hubble |residual] 
(Fig. El. One might well ask whether these SNe, hosted 
by relatively ancient galaxies with little remaining dust, 
also require a nonzero color coefficient TZ. 

From the data considered in this paper, unfortunately, 
the answer is ambiguous. If SN19930 continues to be 
omitted from the Hubble fits, the remaining 12 SNe that 
are hosted by E-I-E/SO-I-SO galaxies show little preference 
for either of the two values of TZ over the other. On the 
other hand, if SN19930 is included, TZ = 2.59 is preferred 
over TZ = 0. Because a robust result cannot hinge on 
including or excluding a single SN, a definite conclusion 
on the color correction of SNe in elliptical, E/SO, and SO 
galaxies is not reached here. 

8. SUMMARY AND DISCUSSION 

Applying both a substantially extended light-curve- 
fitting algorithm and the CMAGIC prescription of Wang03 
to published B and V fluxes from a systematically cho¬ 
sen sample of nearby SNe, we have compared the out¬ 
puts Umax and Bbv of these two types of fit both to 
each other and to Hubble’s-law expectation. After ap¬ 
plying standard corrections for observed SN color, but 
before applying decline-rate corrections, the difference 
Hniax — Bbv reduces to TZAE, where 77. is a generalized 
lin ear co l or-cor rection coefficient of the type measured 
bv lTripnl l|1998l) : and AE = £ — (i?max—Umax) is the dif¬ 
ference between standard measures oi B — V color used, 
respectively, by the CMAGIC and light-curve fit methods. 

After AE is corrected for SN decline rate Amis, close 
agreement in the two measures oi B — V color is achieved: 
their mutual rms is only 0.053 ±0.006 mag. The required 
correction is compatible with a bilinear but not a linear 
Amis dependence (Fig. El; the nonlinear term is highly 
significant (5.8(t). Also, to bring the ^ for the fit to 
AE = 0 down to an acceptable level, the propagated 
photometric errors in AE must be augmented in quadra¬ 
ture by an intrinsic noise term uab = 0.029 ±0.007 mag 
rms. 

In CMAGIC plots of B vs. B — V shown in their Figs. (I- 
4), Wang03 distinguished between SNe that did or did 
not show a “bump” near maximum light. As is evident 
from those plots, SNe which did exhibit a bump had a B 
flux at maximum that was larger than would have been 
obtained by extrapolating to the same epoch the linear 
fit to B vs. B — V that was made in the post-maximum 
CMAGIC region. Therefore the prominence of a bump 
in the CMAGIC plot is intimately related to the differ¬ 
ence between Umax and Bbv, and therefore to the size 
of AE. WangOS noted further that SNe which exhib¬ 
ited a bump tended to be slow decliners, with low Ami 5 . 
The continuous dependence of AE upon Ami 5 depicted 
in our Fig. El suggests that SNe may not be bifurcated 
into bumpless and bumpful subsets; rather, their sys¬ 
tematic differences Umax — Bbv may be described by a 
continuous (but nonlinear) function of SN decline rate 
that reaches its nadir for the slowest decliners, as would 
be expected from Wang03’s observations. 

Turning to results from our light-curve fits, we have 
studied the group of SNe which inhabit the “blue ridge” 
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region |-Bmax — KnaxI ^0.1. In particular, on a plot 
of Smax —Imax VS. Ami 5 (Fig- dl, we souglit but did 
not find evidence for a linear Amis-dependence of the 
B — V color about which this group is clustered: within 
an uncertainty of «0.06 the best fit Amis coefficient A 
is consistent with zero. Taking ^ = 0 as a fair repre¬ 
sentation of the blue ridge topology on this plot, and 
using the result of our best fit to the Amis-dependence 
of AE, we can deduce the topology of the blue ridge on 
a plot of CMAGIC color £ vs. Amis. On such a plot, for 
Amis < 1-1 the color of the ridgeline would redden in 
proportion to Amis with a slope aAE — cx'AE — h-53 (Ta¬ 
bled; at higher decline rates it would deredden with a 
much gentler slope aAE + ol'ae ~ ~0-07. On their plot 
of £ vs. Amis, Wang03 fit the blue ridge to a linear de¬ 
pendence on Amis over its full range, obtaining a slope 
equal to 0.25 ±0.04, in agreement with our average slope 
o^AE = 0.23 ± 0.03. To sum up, when measured at peak 
luminosity the blue ridge color shows little dependence 
on decline rate; but when the color assigned to SNe in the 
ridge is strongly influenced by post-maximum values of 
B — V, as occurs for the CMAGIC color £, the ridge color 
tends to redden as the decline rate steepens, particularly 
for slow decliners. 

From their study of Umax —Imax vs. Amis, Phil99 re¬ 
ported a blue ridge slope A = 0.114 ± 0.037 for SNe 
that were preselected to exhibit low B — V color excess 
in the nebular phase (30-90 days after maximum light). 
This Amis dependence is intermediate between the de¬ 
pendences we observe, with no comparable preselection, 
for Umax —Vmax and for £. Our Amis dependence for 
-Bmax ~ Ftnax is not significantly different from that of 
Phil99. 

A Hubble fit to the average of Smax and Bbv, when 
combined with the fit to AE = 0, is equivalent to a 
global fit that specifies the color and decline-rate cor¬ 
rections for both Hmax and Bbv- It yields the B — V 
color-correction coefficient 77 = 2.59 ± 0.24. This value, 
which parametrizes the effects both of extinction by host- 
galactic dust and of intrinsic SN color variation, lies 
well below the coefficient Rb ~ 4.1 that convention¬ 
ally is used for Milky Way dust. The fit result would 
increase by A77/77 = ±0.054 ± 0.027 if we were to un¬ 
fold the effects of finite color resolution, and it would 
diminish by A77/77 = —0.102 if we were to bias it 
by performing a stepwise rather than a continuous 
minimization, as is discussed in dEHl Our measure¬ 
ment of 77 is consistent with the values 2.09 ± 0.38; 
2.4 4; 2.46 ± 0.46 ; and 2.19 ± 0.33 obtaine d, respectively, 
bv iTripul (1199811: iTriop fc Branchl (1199911 : IParodi et al.l 
(EiOfll : and lOuy ct al .1 (1200,511. (No t e, ho wever, that 
iTripp fc Branchl (1199911 and iGuv et'H] (1200511 report hav¬ 
ing performed a stepwise minimization.) Our result 
for 77 is incomp atible with t h e valu es 1.7 1 and 3.5 used, 
respec tively, by ISaha et HI (1199911 and lAltavilla et al.1 


That 77 is significantly below the canonical value of 
« 4.1 for Milky Way dust is a re al puzz l e that we would 
like to understand. Very recently. IWand (1200, 5ll advanced 
an explanation of the low value oi Rb that involves scat¬ 
tering by dust clouds located in the circumstellar (CS) 
environment. From CS dust at the radius (10^® cm) mod¬ 
eled in that paper, scattered light would be delayed sig¬ 
nificantly relative to unscattered light emitted directly 


by the SN, distorting the evolution of its observed lumi¬ 
nosity and color. In the present paper, where statistical 
analysis of our 56-SN sample yields an example of the low 
value oi Rb which CS dust purports to explain, similar 
statistical techniques may be used to search for such dis¬ 
tortions. For example, larger amounts of CS dust could 
make the decline rate parameter Amis appear smaller as 
a result of the delay in scattered light, causing a negative 
correlation between color and decline rate. 

Pursuing this point. Fig. 3(a-b) of (lWanell200,5ll illus¬ 
trates the effects of that paper’s CS dust model on a 
single SN. Its B — V color is 7? « —0.03 when unextin¬ 
guished, vs. E « 0.24 after absorption and scattering 
by CS dust. Taking into account the change in time 
of 7?i„ax, Amis drops from « 0.95 without extinction 
to « 0.70 after CS dust extinction, yielding the quo¬ 
tient Q = A{Ami^)/AE « -0.25^.27 « -0.92. If 
one interprets Q as the coefficient of a linear depen¬ 
dence upon E of (observed) Amis, for these variables 
one predicts a Pearson correlation coefficient Rpred = 
Q a{E)/a{Ami^) ~ —0.52 ±0.08, where the error arises 
from uncertainties on the standard deviations a oi E 
and Amis. For comparison we examined the 43 late- 
type SNe shown by closed triangles and closed circles 
in Fig. 12 whose correlation between E and Amis is 
R-obs = 0.07 ± 0.16. Barring a > 3cr fluctuation, Robs 
would have been sensitive statistically to the CS dust 
prediction. However, if part of the spread in E were due 
to variation in intrinsic color Aintr, and if true Amis were 
positively correlated with Aintr (as suggested by Fig. 5(a) 
of Phil99), a negative correlation due to CS dust could 
have been cancelled, accounting for our low value of Robs- 
We also searched for a more general dependence of Amis 
upon E by dividing the same 43-SN sample into reddened 
vs. unreddened subsamples (filled circles vs. filled trian¬ 
gles in Fig. 12. Their di fference in (E ) matched AE of the 
single SN modeled by (IWangll20?)5ll . Again, at a compa¬ 
rable level of sensitivity, we did not observe the relative 
drop in (Amis) predicted for the reddened subsample. 
As models of circumstellar dust around SNe are refined, 
statistical constraints of the type described here should 
prove valuable. 

Resuming the summary of this paper’s results, we re¬ 
turn briefly to SNe that inhabit the blue ridge. There we 
selected a subsample of SNe that offer a relatively good 
chance of not requiring any color correction. Notwith¬ 
standing the possible bias in their selection, in separate 
fits to the Hubble residuals of 7?max these SNe prefer 
77 = 2.59 to 77 = 0 with a likelihood ratio of 1000:1. 

The nonlinear decline-rate correction required to rec¬ 
oncile the 7?max —Vmax and CMAGIC colors forces Hmax — 
Bbv^ which is proportional to the color difference, to ex¬ 
hibit the same nonlinearity. Therefore, with the same 
5.8 ct significance, a nonlinear decline-rate dependence 
must characterize either 7?max, Bbv, or both. When 
the global fit is used to apportion that nonlinearity be¬ 
tween B^ax and Bbv, the nonlinear term is insignifi¬ 
cant in Bbv (l-2cr) and is only marginally significant 
in ilmax (2.7cr). Therefore we conclude only that the 
required nonlinearity is unlikely to be confined fully to 
Bbv- Accepting our best-fit nonlinear parametrization 
of the Amis-dependence of 7?max (Fig.©, over the indi¬ 
cated range of Amis its maximum |deviation| from the 
best-fit linear parametrization is 0.13 mag. Therefore 
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the need for a nonlinear decline-rate correction should 
not seriously challenge the validity of cosmological re¬ 
sults based on past high-z SN observations, upon which 
scant photon statistics have imposed other uncertainties 
that typically are greater. 

Phil99 introduced a decline-rate correction that was 
quadratic in Amis. In Hubble fits to the more limited 
data then available, their quadratic term was not signifi¬ 
cant (0.9cr), and the sign of their best fit nonlinearity was 
such that faster decliners require steeper corrections. In 
the present analysis, slower decliners require steeper cor¬ 
rections (Fig.CJ. However, it is unreasonable to compare 
Fig.ddirectly with the results of Phil99; their Hmax mag¬ 
nitudes were corrected for host-galactic reddening using 
a color coefficient Rb « 4.1 that is characteristic of Milky 
Way dust, whereas ours were corrected for host-galactic 
reddening plus intrinsic SN color variation using the co¬ 
efficient TZ = 2.59 ± 0.24 that we measure in the current 
analysis. We note that, for slow decliners, our best fit 
Ami 5 correction to the CMAGIC output Bbv is consis¬ 
tent with nil. 

The proportionality of Bmax — Bbv to the color dif¬ 
ference /S.E causes the measured noise ttab to map into 
a mutual intrinsic scatter uab = 0.074 ± 0.019 mag in 
Hmax vs. Bbv- The values of for Hubble fits to Hmax 
and Bbv are such that the maximum likelihood occurs 
when all of ttab is confined to Bbv- However, the data 
do not significantly prefer that solution to one in which, 
for example, 1/^/2 of uab is apportioned to each lumi¬ 
nosity measure. At the level of a 44:1 likelihood ratio, it 
is unlikely that all of bab is confined to Hmax- Even the 
full 0.074 mag is not a high level of scatter compared to 
the additional errors of « 0.15-0.18 mag often added in 
quadrature to photometric errors on the luminosities of 
high-z SNe. 

Of course, additional intrinsic scatter that similarly af¬ 
fects both Hmax and Bbv may be present. Lacking the 
statistical power to separate it from the effects of pecu¬ 


liar velocity, we lump those noise sources into an effec¬ 
tive rms SN velocity z;eff = 382^*^52 km s“^. This places 
a 95% confidence upper limit of 486 km s“^ on the true 
rms line-of-sight peculiar velocity of nearby SNe relative 
to the Hubble flow. Most recent SN analyses employ 
values of the rms peculiar velocity that are below this 
limit; a particular except i on is t he recent SN calibration 
paper of lWang. X. et all l|200,5ll . whose analysis used a 
value of 600 km s“^. An rms peculiar veloci ty of 500 km 
s~^. b arely above our limit, wa s assumed bv lTonrv et alJ 
(l2f)(m an d IBarris et all (1200411 . 

Applying three different measures of dispersion, we 
find that Bmax, Bbv, and their unweighted average 
are scattered about the Hubble line by amounts that 
are comparable and of order 0.14 ± 0.02 mag (for 
0.015 < z < 0.1), and that are not easily distinguished 
from each other (Table [TJ. The uncertainties inherent 
in such measures are substantial, not only from random 
fluctuations but also from systematic variations in iden¬ 
tification and treatment of outliers and in other method¬ 
ologies. We urge caution in interpreting these and other 
estimates of the uniformity of Type la SNe. 

Our priorities for extension of the work described in 
this paper include achieving a better understanding of 
correlations among the published uncertainties in raw 
flux measurements made on individual SNe, and putting 
that understanding to use in further-refined light-curve 
fits. In particular, we look forward to better quantifica¬ 
tion of the systematic errors that are assigned to the fit 
outputs. 


We thank Weidong Li for providing photometric data 
for SN2002el prior to publication. This work was 
supported by the Director, Office of Science, of the 
U.S. Department of Energy under Contract No. DE- 
AC02-05CH11231. 


APPENDIX 


A. STATISTICS 

In the Introduction, reference is made to the possibility of a general solution to the problem of estimating the 
deviation A of a single SN from a fiducial absolute SN magnitude, using a set x = {xj} {l<j<n) of its measured 
parameters. Given a training set consisting of an ensemble of N SNe in the nearby Hubble flow whose corresponding 
deviations A^ (l<i<N) are deduced from their redshifts, an example of such a general solution is provided by the 
(JPDE recipe of lKniiteson et ahl (1200211 . A covariance matrix 


N 




(Al) 


is constructed, where = (x^, A^), is the mean v^, and k,l index the (b-|- 1) components of v. The best estimate 
of A maximizes the joint probability density 


p(x, A) 


I 

iV(v/^/i)("+i)deti/2(S) 


N 

Xexp 


(v-vQlS ^(v-vz) 
2h^ 


(A2) 


where h is a smoothing parameter determined by N and n. 

In H6.I1I reference is made to the unweighted rms Hubble residual a. It is defined by 


(7^ = 


N 


Ndof ^ 
1 — 1 


(rrii - 


(A3) 
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Day after B maximum 

Fig. B1.— Six-parameter light-curve fits to B and V vs. rest-frame epoch for the highly extinguished Type la SN 1995bd. Measured 
data points are shown by circles, with error bars indicating the Icr uncertainties. 


where A^dof is the number N of SNe with corrected magnitude rui at redshift Zj, less the number of parameters in the 
ht function Reference is made also to the weighted rms (wrms) Hnbble residual a-w. It is defined by 




N 

fVdof 


N 

E 

2=1 


{mi - fi{zi)f 


N 

'E: 

2=1 


(A4) 


where Ui is the total uncertainty in taking into account photometry, intrinsic noise, and peculiar velocity. 


B. LIGHT CURVE FIT DETAILS 

In H2.ll reference is made to the “stretched” SN phase t* introduced in equation d). We used the function 

t (t;so,ro)=t(^l + (—• (Bl) 

Also in 30 reference is made to the fact that our light-curve hts to a minority of SNe employed six parameters 
(including the delayed parameters eo and td in equation d)) rather than four. As is seen in Table El 18 of 61 htted 
SNe fell into that category. Of those, 16 are the subset of the 61 SNe that were identified spectroscopically to be 
SN1991T-like. As well, adding two parameters improved the quality of ht to SN1996X and SN1998bu. Figure IHTI 
exhibits a sample 6-parameter light-curve ht. 

In dif is asserted that the quantity E (equation d) roughly tracks the photometric error on AE = £ — E. From 
the dehnitions of £ and E in m 

PbvAE = PbvV^Z - i^BV - . (B2) 

After approximating I3bv ~2, neglecting uncertainties in corrections for Milky Way extinction, and taking the three 
raw quantities in equation (IB2II to be uncorrelated (see H2.5II . the photometric error on its right-hand side reduces to 

E. 

Also in m reference is made to adjustment of the photometric ht errors that are hagged by the footnote in Table 
El These adjustments were occasioned by the detailed graphical inspections of ht results noted in 33 When the result 
of a photometric ht showed a systematic deviation from measured huxes, usually owing to rigidity of the light-curve 
template, we left all best-ht values untouched but perturbed upward the error assigned to the ht parameter, e.g. 
to bring it into < 2cr agreement with the data. No account was taken of the implications for Hubble hts. When the 
error in B^'^^ was scaled, so was the error in Amis. Overall, 13% of the errors quoted in Table El were adjusted. 

C. EXAMPLE OF Hmax = Bbv 

Under the special conditions enumerated in 3231 Hniax and Bbv are fully correlated. As a simple example that 
satishes these conditions, suppose that only three measurements {Bi,B 2 ,B-i] of raw B magnitude are made. In the 
LC ht to the B band, after one hxes tq = oo and eg = 0 in equation o and optimizes the three remaining parameters 
Fq, to: and sg, normally one obtains a best ht light curve that passes exactly through the three measured B points. 
Suppose further that only a single measurement Vg of U magnitude is made, at the same epoch tg as Hg. In the LC ht 
to the V band, to and Sg must then be determined by their H-band values, while, in this limiting case, only Fg may 
be optimized. Again the best ht light curve passes through the single measured V point. Within an additive constant, 
using the fact that all residuals vanish, 

Bi::^ = Bo-mB{n)- (Cl) 

V^Z = V3-mv{n); (C2) 

E = Bg-Ug-(mB(rg)-mv(Tg)) , (C3) 
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TABLE D1 

K CORRECTIONS TO Bmax, Vmax, Amis, BbV 

(mag), and Pbv, for four SNe with redshifts 

AND COLORS LISTED IN THE LAST TWO ROWS. 


K correction 


for SN 
to quantity 

1992bp 

1992al 

1995ak 

2000cn 

-Smax 

-0.027 

0.009 

0.008 

-0.013 

Vmax 

0.043 

0.012 

0.009 

-0.005 

Amis 

-0.118 

-0.019 

-0.049 

-0.093 

Bbv 

0.223 

0.022 

0.037 

0.040 

Pbv 

-0.226 

-0.019 

-0.038 

-0.038 

Characteristics of iC-corrected SNe: 

z 

0.0786 

0.0135 

0.0230 

0.0232 

E 

-0.050 

-0.051 

0.000 

0.202 


TABLE D2 

Coefficients of 2 , 2 ^, an d 2 ® from the cubic 

FITS SHOWN IN FlC. IDlf fei. (c), (d), AND (e). 


Fig- 

Im1 

Coefficient 
of 2 

Coefficient 
of 2 ^ 

Coefficient 
of 2 ^ 

(6) 

-0.97119 

28.9034 

-187.314 

(c) 

3.86474 

-23.3410 

165.731 

(d) 

-1.17487 

16.5470 

-259.775 

V) 

3.04291 

-10.4068 

149.714 


where ms.v are the template magnitudes and T 3 = (t^ — to)/so- If ^3 lies within the CMAGIC linear region, one may 
fix I3bv at a representative value /3o and use the single pair {Ra, V 3 } to determine ; within a constant, 

i?W=/3oI"3-(/3o-l)S3 ; 

Forming AE = £ — E and taking its differential, 

d(A£') = d(mB(T 3 ) - mvlTs)) - dmBiT3) . (C 6 ) 

If the templates are “magic” - that is, if they satisfy the condition diriB/d{mB—iTLv) = Po in the CMAGIC linear region 
- the differential d(Ai?) vanishes: apart from an additive constant, the two color measures E and £ are equivalent, 
and so are Rmax and Bbv (eqs. 0 and ©)• Therefore, if the data are minimal and the templates are magic, the 
maximum-luminosity and CMAGIC analyses are mutually equivalent. 

D. K CORRECTION DETAILS 

A brief description of our approach to K corrections is found in 01 Table inn exhibits the difference between 
A-corrected and A-uncorrected Umax, Emax, Amis, Bbv, and Pbv, for four selected SNe whose redshifts and colors 
E = B^^^ — are provided in the final two rows. SN1992bp was selected because it has the largest redshift in 
the sample; SN1992al differs from it primarily in redshift. SNI995ak and SN2000cn were included because they differ 

from each other primarily in color. _ _ _ 

Additional details of the A correction method are provided in Fig. IdTI and Table El Figure IrTiT a) exhibits linear 
fits to the dependence upon B — V color of the A corrections to B and E magnitudes of SN1994D. The Icr error bars 
are estimated by varying the choice of library spectral template while warping it to reproduce the observed SN color. 
The slopes and inter cepts of these lines are fit to cubic functions of z shown in Figs. IPlT b-e). whose fit coefficients 
are provided in Table El The A corrections to Rmax, Enax, Amis, Bbv, and Pbv are obtained by evaluating the A 
corrections to B and/or E magnitudes at the appropriate B — V colors. 


(C4) 

(C5) 


E. NOISE DETAILS 


In gaussian approximation and in the basis {i?max, Bbv} of the measured B magnitudes, the noise in a single SN’s 
measured luminosity from all sources other than photometric error is expressed by the noise covariance matrix 


Y = 


max 

Ee 



(El) 


where the off-diagonal element E in principle is bounded only by —o’maxO'Bv < E A CTmaxCsy. 
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obs obs 




Z Z 




Fig. D1.— Illustration of K correction method, (a) K corrections at ^ = 0.05 to B magnitudes (full circles) and V magnitudes (open 
circles) of SN1994D vs. observed B—V color. The lines are linear fits with intercept fi and slope i/i. (b-e) Dependence on z of {b,d) fi and 
(c,e) for (fe,c) B and (d,e) V magnitudes. The fit curves are cubic polynomials in z. 


Because allocating the measured noise between B 
approximate Y as the sum of two simpler matrices 


Y = 


max and Bbv is challenging statistically, it is advantageous to 


) 

^u,BV J 




(E2) 


The first covariance matrix arises from sources of noise cr^ that are completely uncorrelated between Bmax and Bbv ^ 
while the second represents noise sources (including peculiar vel ocity and, at higher redshift than con sider ed here, 
weak leasing) that affect both measurements equally. Equation is an approximation to equation m only in 
that it restricts the off-diagon al ele ment of Y to the range 0 < 14 < min (Umax) ^%v)- (This is the “ mild assumption” 
to which reference is made in lO I Intrinsic noise c an co ntribute to either of the terms in equation 1E2|I . 

After the noise covariance matrix Y from equation (IE2II is transformed to the basis {AS, (B)}, its diagonal elements 
are 


2 2 2 

= ^u.max + ^u,BV '■> 

^\b) = K'^u.max + ^u,Bv) + ■ 


(E3) 

(E4) 


Apart from sources of noise 14 that are common to both Bmax and Bbv, it is evident from eqs. (IE:IIE4II that the 
rms noise uab that is measured in the comparison of Bmax to Bbv propagates into (B) as n(^B) = n-Ae/S. 
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In >15.21 reference is made also to the possibility of using the redshift dependence of the common-mode noise to 
resolve the effects of peculiar velocity from those of other sources. Such analysis, which is not pursued here due to 
insufficient statistics, would use an equation of the form 


^2 _ 2 

(T — (T 

max u,max 




int 


C V, 
z 


pec 

C 


(E5) 


where is the portion of 14 that is ^-independent. The factor C{z) is equal to 5, multiplied by a correction that, 
for a flat universe with ^1m = 0.3, ranges from 1.00 to 1.05 over the redshift range (0.003 to 0.079) of the actual SNe 
in our fitted sample. 


F. COLOR RESOLUTION DETAILS 

In >15.51 reference is made to two different methods of estimating the correction to TZ from the effects of finite color 
resolution. In these methods the average color (E) (equation (I12II 1 was used. In the first method, the SNe were 
redistributed in color so that a twice-resolution-smeared color ideogram of the redistributed SNe matched a once- 
resolution-smeared ideogram of the unredistributed SNe. During the redistribution, the SN excursions were limited by 
requiring x^/dof = I, i.e. the rms excursion of a SN was limited to one standard deviation in its own color resolution. 
Using the color-redistributed SNe as indicators of the true SN color distribution, a linear least-squares fit was made 
to a scatter plot of redistributed color vs. observed color for each SN. The fit slope was 0.962, implying that when TZ 
is fit to the observed SN color it is underestimated by 3.8%. 

In the second method, a once-resolution-smeared color ideogram first was mildly smoothed using a multigaussian 
to yield an “apparent” color distribution. Next, a continuous “true” color distribution was identified such that after 
two gaussian smearings it matched the apparent distribution. In this method the smearing width was set to the 
median color resolution of the SN set. Next, by smearing the true distribution only once, the average true color that 
contributes to each apparent color was determined. Then an analytic function was fitted to a plot of (average true 
color — apparent color) vs. apparent color. Using this function, the ratio p of average true color to apparent color 
was constructed as a function of apparent color. Finally, the global fit described in was repeated with each SN’s 
observed color multiplied by p. Relative to the standard ht value, TZ was found to increase by 7.0%. 
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